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Abstract 

The 2.5(7 discrepancy between theory and experiment observed in the difference 
Acp(B~ — > ir°K~) — Acp(B° — > ir + K~) can be explained by a new electroweak pen- 
guin amplitude. Motivated by this result, we analyse the purely isospin-violating decays 
B s — > <pp a and B s — > 07r , which are dominated by electroweak penguins, and show that 
in presence of a new electroweak penguin amplitude their branching ratio can be enhanced 
by up to an order of magnitude, without violating any constraints from other hadronic 
B decays. This makes them very interesting modes for LHCb and future B factories. 
We perform both a model-independent analysis and a study within realistic New Physics 
models such as a modified-Z°-penguin scenario, a model with an additional Z' boson and 
the MSSM. In the latter cases the new amplitude can be correlated with other flavour 
phenomena, such as semileptonic B decays and B s -B s mixing, which impose stringent con- 
straints on the enhancement of the two B s decays. In particular we find that, contrary 
to claims in the literature, electroweak penguins in the MSSM can reduce the discrepancy 
in the B — > TtK modes only marginally. As byproducts we update the SM predictions 
to Bt{B s -> 07T°) = 1.6±o 3 • 10" 7 and Br(B s -> <f>p°) = 4.4±§; 7 . • 10~ 7 and perform a 
state-of-the-art analysis of B — > wK amplitudes in QCD factorisation. 
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1 Introduction 



At present flavour physics has entered a new exciting era. The new experiment LHCb and the 
planned super-B-factories will bring the precision of Standard Model (SM) tests and the scope 
of searches for New Physics (NP) to unseen heights. Particularly important thereby are flavour- 
changing neutral current (FCNC) decays, which in the SM are highly-suppressed electroweak 
loop processes. In this work we present a phenomenological analysis of two hadronic FCNC 
decays, namely B s — > (ftp and B s — > (fttr . We argue that within the next years these decays will 
become very interesting objects for experimental analyses of the electroweak penguin sector. 
Up to now, this sector has been tested in hadronic decays only in B — > ttK modes, and the 
discrepancies found between the SM prediction and experimental measurements is the main 
motivation for our work. 

The four B — > ttK decay channels, first observed by the CLEO experiment in the late 1990s 
0, Q], have become by now a classic in flavour physics thanks to the precise measurements by 
BABAR and BELLE. This is also reflected in the large number of theoretical studies of these 
decays in the SM and various extensions of it. Charged and neutral B mesons can decay to 
a ttK final state due to a weak process at the partonic level, b — > sqq with q = u,d. This 
process is dominated by an FCNC loop governed by the CKM factor V t *V t b and receives, in the 
q = u case, also a small tree-level contribution involving the smaller CKM factor V* s V u b- The 
B —7- nK branching fractions are therefore small, of order 0(1O -6 ), and sensitive to new FCNCs 
arising in extensions of the SM. For this reason they are, together with the corresponding CP 
asymmetries, important observables for tests of the SM flavour structure and for NP searches. 

With the data of the B factories having become more and more precise, some discrepancies 
between B — > ttK measurements and SM predictions have occurred, provoking speculations on 
a U B — > ttK puzzle". To date, the measurements of the branching fractions have fluctuated 
towards the SM predictions, the latter still suffering from large hadronic uncertainties, and 
only the CP asymmetries show an unexpected behavior manifesting itself in the quantity 

AA CP = A C p(B- it K-) - A C p(B° ir + K~). (1) 

For this observable we find in the framework of QCD factorisation (QCDF) 

AA C p^ 1.9+11% (2) 

as the SM prediction, which differs significantly from the experimental value 

Ai C p= P ' (14.8 ±2.8),%. (3) 

Adopting a frequentist approach where we consider a theoretical "error bar" as a range of 
values definitely containing the true theory result but without assigning any statistical meaning 
to it 0], this amounts to a 2.5a discrepancy. 



A point which has received much attention in the literature (see e.g. |8[ and references 
therein) is the fact that the formerly observed discrepancies as well as the currently existing 
anomaly in A^lcp suggest a violation of the strong isospin symmetry beyond the amount 
expected in the SM. This has often been interpreted as a hint for enhanced electroweak penguins 
(EW penguins) 0-11]. We will give a brief overview and discuss the current status of this topic 
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in section I2TT1 Whether the 2.5a discrepancy in A^4cp is a hint for NP in EW penguins or a non- 
perturbative hadronic effect or simply a statistical fluctuation is controversial. The point that 
we want to make is that, in order to assess this question, it is highly desirable to obtain further 
information from other hadronic decays which are sensitive to EW penguin contributions. For 
this reason we study the purely isospin- violating decays B s — > <pp° and B s — > (pir , which are 
dominated by EW penguins, extending and updating our analysis presented in ref. If NP 
in this sector exists at a level where it can explain the A^4cp puzzle, it could be clearly visible in 
these purely isospin-violating decays. The upcoming new generation of flavour experiments will 
have the opportunity to detect these modes for the first time and to measure their branching 
fractions. The aim of our work is to provide a detailed analysis from the theory side, both in 
the SM and beyond. 

Since the decays B s — > 4>p° and B s — > cp-K are not related to other decay modes via flavour 
symmetries, the non-perturbative part of their decay amplitudes has to be determined from 



first principles. This can be achieved using the framework of QCDF 13l4l7l|. This method 
amounts to a calculation of the hadronic matrix elements up to corrections of order Aqcd/^B; 
where Aqcd ~ 0(200 MeV) is a typical non-perturbative energy scale of strong interactions. 
We will use this method throughout the paper in all analyses of B decays to light mesons. 

The plan of the paper is as follows: In Chapter [2j we discuss the issue of isospin- violation 
in B — > ttK decays and the phenomenology of B s — > (j)p° and B s — > (pn . As a byproduct we 
provide simple formulas which allow for an easy calculation of various observables concerning 
these decay modes, taking into account NP effects in EW penguins. Chapter [3] contains a 
detailed quantitative analysis of B s — > (j)p and B s — > 07r° in different scenarios of a model- 
independent parameterisation of NP in EW penguins. This analysis is performed in light of 
our present knowledge on EW penguins from other B decays, in particular B — > irK. It 
is complemented in Chapter [J] with studies of particular extensions of the SM which feature 
enhanced EW penguins. We conclude in Chapter [5l We keep the main body of the paper free 
of technicalities and refer the reader interested in technical details to the appendices. 



2 Isospin-violation in hadronic B decays 
2.1 The B ->- tt/\ modes 

The B — > ttK decays are dominated by the isospin-conserving QCD penguin amplitude. Nev- 
ertheless, they receive small contributions from the tree and the EW penguin amplitude, which 
are isospin-violating. Combining measurements of the four different decay modes B~ — > 7r~K°, 
B~ — > ir°K~, B° — > ir + K~ and B° — > ir°K°, it is possible to construct observables in which 
the leading contribution from the QCD penguin drops out, so that they are sensitive to isospin 
violation. 

The mesons participating in B — > irK decays transform under isospin rotations as 

(B°,-B-) 1/2 , (K°,-K-) 1/2 , (fr+.-^-Tr-)!. (4) 
Furthermore we can assign isospin to the operators appearing in the effective Hamiltonian 

HcS = 7f E ^ ( C ^l + C ^ + E C ^ + C ^7 7 + C 8g Q 8 g\ + h.C, (5) 
v p=u,c \ i=3 / 
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(s) 

which mediates the B — > nK transitions. Here X p = V p bV* s represents a product of elements 
of the quark mixing (CKM) matrix, Q\ 2 are the so-called current-current operators, Qz,...,& 



are QCD penguin operators, Qj^ and Q$, 
operators and 



represent the electromagnetic and chromomagnetic 



Q 7 = (s a b a ) V - A J2 q %eq(qpqp)v+A, 
Q 9 = (s a b a ) V - A J2 q \eq{qpqp)v-A, 



Qs 

Qio 



W-A 



E 



q 2 e <j( 



iqa)V+A, 



(6) 



(s a b ) V -AEn %eq(qpQa)v-A, (7) 



are the EW penguin operators (a, j3 denote colours). The latter are of great importance for 
our work. We define the operators as in 15(] so that Ci(M\y) = 1 at leading order. Containing 
uu- and dd-bilinears, the operators Qi,...,Qio can be distributed among 



H 



eff 



n,AI=l | n,AI=0 

rt eff ' n cS 



(8) 



according to the decomposition 1/2(8)1/2 = 1 © [la ]. Since the QCD penguin operators 
Qi,.. .,6 involve the isosinglet combination (uu + dd), they contribute solely to T-L^ =0 whereas 
the other operators give contributions to both parts of % e ff- The B — > ttK decays thus follow 
the isospin pattern 

A T — i n 

3/2 1/2, (9) 



1/2 A i=l'° 1 ® 1/2 



implying that all four decay amplitudes can be decomposed into three independent isospin 
amplitudes, -4f/2 = °> "^i/2 =1 anc ^ -^3/2 =1 with the lower index denoting the total isospin of the 
final state. 

One finds that B — > ttK is dominated by the QCD penguin contribution and thus |^.^2 = °l ^ 
l"^3/2 = i/2l' ^° a ^ rs ^ a PP rox i ma ti° n ) all the decay modes can be described by the amplitude 
*^i/2 =0 on ^' dictating the relative size of the branching fractions to be 1 : 2 : 1 : 2 (in the same 
order as in tab. [I]). 

The isospin-invariant amplitudes receive contributions from various SM quark diagrams. 
It is only at the level of these diagrams that the pattern of CP violation can be correctly 
implemented, i.e. that the amplitudes Ay 1 ^ , *^3/2 = i/2 can be related to their CP-conjugated 

counterparts Aij 2 °, ^3/2, 1/2 • This suggests an alternative parameterisation of the amplitudes 
in terms of the topologies of the underlying quark-level transitions [191 . 1201 ]: 



A(B~ - 
V2A(B- - 
A(B°- 
V2A{B° 



1 



tt-K") ~ P I 1 - ^r£ w + =-ri w 



tt u K ) ~ P ( 1 + r EW + ^ew + ^r£ w - (r T + r c )e 



-17 



TT+K-) ~ P 



2 c 

1 + g r EW 



PI- r E w 



- (7 



(10) 



This topological parameterisation is illustrated by the corresponding Feynman diagrams for 
B~ — > -k®K~ in fig. [TJ In eq. (fTUj) we have factored out the dominant QCD penguin amplitude 
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b W s u b W s 




Figure 1: Diagrams representing the topological parameterisation in eq. (fTOj) for B~ — > ir°K~. First 
line from left to right: QCD penguin (P), colour- allowed EW penguin (rgw), colour-suppressed EW 
penguin (r£ w ). Second line from left to right: colour-allowed tree (rx), colour-suppressed tree (re), 
EW penguin annihilation (?"g W ). 



P and neglected penguin amplitudes suppressed by IKTs^&I/IKs^c&I- The dependence on the 
weak CKM phase 7 has been made explicit, while strong phases are contained in the ratios 
Ti which fulfill I I < 1. These quantities denote corrections from different types of Feynman 
diagrams: rx and tq stem from colour-allowed and colour-suppressed tree diagrams, tew and 
r^ w from colour-allowed and colour-suppressed EW penguins, respectively. Annihilation via 
QCD penguin diagrams is absorbed into P whereas weak annihilation via EW penguin diagrams 
is parameterised by r^ w and colour-suppressed tree annihilation is neglected. With our QCDF 
setup explained in Appendix [A] and the expressions for the ratios given in Appendix [B] we 
obtain 

rrr _ 1 7 +0.07 , o no+0.03 • 
'T — u - x '-0.06 ' u - uo -0.10*' 



rc = 0.07±85S + ("0.01)1831 



r E w= 0.13^ + 0.02^ 2 7 z, 
rgw = 0.04i^ + (-0.01)1^ i, 

r^ w = 0.007iHSo + (-0.004)1°;^^. (11) 

The result displays the typical features of QCDF predictions, namely small strong phases 
and large uncertainties of colour-suppressed topologies. The smallness of the r% reflects the 
domination of the isospin-conserving QCD penguin and justifies the expansion of physical 
observables in the rj. Among the isospin- violating contributions the colour-allowed tree gives 
the largest corrections followed by the EW penguin which dominates over the colour-suppressed 
tree. The colour-suppressed EW penguin ratio r^^y and especially the EW penguin annihilation 
ratio r^ w are quite small and consequently they have been omitted in most analyses of B — > ttK 
decays. In particular, the possibility of having NP in the EW penguin annihilation amplitude 
r^ w has to our knowledge not been considered so far. However, we want to point out that 
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Table 1: Theoretical vs. experimental results for B 
Q • The original results can be found in 21 - 35| . 



ttK decays. The experimental data is taken from 



such an approximation is not valid in the analysis of CP asymmetries: non- vanishing direct CP 
asymmetries are caused by the interference of parts of the decay amplitude with different weak 
and strong phases. Consequently direct CP asymmetries in B — > ttK cannot be generated by 
the QCD penguin amplitude alone and are automatically sensitive to subleading contributions, 
encoded in the imaginary parts of the coefficients. These, in turn, are generated in QCDF 
either perturbatively at 0(a s ) or non-perturbatively at ©(AqcdA^b)- At 0(a s ) the colour- 
suppression of r^ w is not present anymore and the Aqcd /wb - suppressed r^ w can compete 
as well. Therefore we keep r^ w and r^ w in our calculation and we will see in later chapters 
that we can indeed have a large NP contribution in these amplitudes. 

One can easily see from eqs. (JlOllip that the two amplitudes involved in A.Aqp differ only 
by the subdominant contributions rc, ^ew an d r^, all of which are isospin-violating. Turning 
to the CP asymmetries, one finds in the SM 

Aqp(B~ — > 7r°A~) ~ — 21m (rx + rc) sinj, 

A C p(B° -> tt + A~) ~ -2Im(r T )sin7, (12) 

with terms quadratic in the r« being neglected. Thus the only possible explanation for a large 
A^4cp i n the SM seems to be a large imaginary part of rc, i-e. a large absolute value and large 
strong phase of the colour-suppressed tree amplitude, generated by some hadronic effects at 



5 



the low scale Aqcd which can hardly be calculated perturbatively. However, QCDF predicts 
only a small Im(rc), insufficient to explain the data, even when all the theory uncertainties 
are included. Therefore one is tempted to conclude that the discrepancy in A^4cp is not due 
to our lack of understanding of strong interactions but due to isospin-violating NP. 

For this reason, A^4cp has been studied in various NP models in recent publications 0, [36r~ 



441 ] . The main ingredient of these analyses is usually an enhancement of the EW penguin 
topologies by effects of virtual heavy particles. Such contributions can be included into the 
amplitudes (fTOj) by the replacements 



rEW^r EW + r" E we rg w -> rg w + fg w e lS , rj| w -> r£ w + r£ w e lS , (13) 

0) 

where 5 is a new weak phase and are complex numbers including a strong phase. The CP 
asymmetries then become 

A C p(B~ -> tt°K~) ~ -2 Irn (r T + r c ) sin 7 + 21m ^f EW + ^fg w + t^ew^ sin *> 
A C p(B° -> 7T+A-) ~ -2Im(r T )sin 7 + 2Im ( -fg w - -f^y ) sinJ, (14) 



3 " vv 3 



such that 



Lcp ^ -21m (rc) sin 7 + 21m (f E w + ? EW ) sin 5 (15) 

can turn out to be much larger than in the SM. The observed discrepancy can be solved by a 
?ew or a fg W term comparable in size to the corresponding SM term rEW- 

Apart from A^4cp one can also construct other observables from the B — > ttK data which 
are sensitive to isospin violation, for example certain ratios of branching fractions. Even though 
tensions with experimental data in these observables raised the formulation of a "B — > ttK 
puzzle" in the first place 0, 0, [ll|, 4^- 4^] , in the meantime these quantities are in reasonable 



agreement with the SM predictions. However, they serve as important constraints for NP in 
EW penguins and we define and discuss them in Appendix [B] Note in particular that the 
quantity AA CP defined there, which is the difference of the two remaining CP asymmetries 
not appearing in Aj4cp, probes the same combination of Im(fEw) an d Im(f^ w ) as AAcp- 
Unfortunately, data on Aqp(B~ — > 7r~K°) and especially on Acp(B° — > n°K°) are not good 
enough yet to gain any information from these observables. Experimental results and SM 
predictions for the B — > ttK observables are given in tab. [T] 

The main problem which makes it difficult to single out a possible NP contribution in 
B — > ttK decays is evident from (|1U|) : the colour-allowed EW penguin contributions and colour- 
suppressed tree contributions enter the amplitudes in ()10p exclusively in the combination 



fEW ~ r c e 11 . (16) 

This implies that colour-allowed EW penguins and colour-suppressed trees are inextricably 
linked with each other, reflecting the fact that the topological parameterisation contains some 
redundancy. Physical effects found in any experiment cannot unambiguously be attributed to 
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one or the other partner of this topology pair. A new EW penguin contribution rEw e %S can 
be probed only in one of the four physical combinations 

Re(fEw) cos <5 — Re(rc)cos7 + Re(rEw) > 

Im(fEw) sin 5 — Im(rc)sin7, 

Im(fEw) cos 5 — Im(rc)cos7 + Im(rEw) , 

Re(fEw)sin<5 — Re(rc)sin7. (17) 

Therefore, probing fEw e ~* 5 is challenged by the large hadronic uncertainties in the QCDF 
prediction for rc, which can mimic or hide such a NP signal. One possible way to constrain rc is 
the approximate SU(3) flavour symmetry which relates it to a corresponding B — > tttt topology. 
Using this symmetry it has been found that current data on CP violation in B° — > tt°Ks is 
also in disagreement with the SM, independently of AAqp, and can be explained by adding 



tew to the amplitude [421 ] 



The perspective of our work is the following: In order to find out whether the AAqp 
discrepancy really is a manifestation of isospin-violating physics beyond the SM, one should 
also study other observables on which such a kind of NP could have a large impact and see 
whether similar effects appear in measurements of these observables. Our proposal in this 
work is to test the hypothesis of isospin-violating NP by looking at processes which are highly 
sensitive to it, namely purely isospin-violating B s decays. 



2.2 Purely isospin-violating B s decays 

EW penguin contributions to hadronic B decays are usually overshadowed by the larger QCD 
penguins. This problem can be avoided if one succeeds in probing exclusively the AI = 1 part 
of the effective Hamiltonian which is orthogonal to the QCD penguin operators. To achieve 
this for B — > irK, we had to single out the AI = 1 part of the transition in eq. © by combining 
different isospin-related decay modes, for example by considering the observable AAqp- Our 
proposal now is to consider decays to which QCD penguins do not contribute at all, i.e. pure 
AI = 1 decays, where no such procedure is needed. 

There are no two-body decays of the B° or B^ meson with this property. In these cases 
the final state would have to be a pure |3/2,±l/2) isospin state which cannot be constructed 
out of two mesons. The B s meson, on the other hand, is an isosinglet and it can decay as 

^ 0(8 1 = 1. (18) 

The final state must consist of an isospin triplet, i.e. 7r° or p°, and an isosinglet, i.e. a meson 
with the flavour structure ss. In order to avoid complications stemming from r] — ry'-mixing, 
we restrict ourselves to the vector-meson (f> which is to a good approximation a pure ss state. 
This leaves us with the two AI = 1 channels 

B s -> (j)p° and B s <f>7r°. 



7 





Figure 2: penguin, tree and annihilation topologies contributing to B s — > <f>p ,4>ir z 



So far only an upper limit Br(l? s — > 4>p°) < 6.17 • 10 -4 exists [49J and no detailed theory 
analysis of B s — > (ftp ,(f)ir has been published. Only the SM branching fractions and CP 
asymmetries have been calculated in general surveys on B decays to light mesons 
In addition, J3o — > 07r° has been suggested as a tool to measure 7 via the mixing-induced CP 



asymmetry [501 ]. Since in the era of LHCb and super B-factories these two processes will become 
interesting objects for tests of isospin-violation and potential NP we will in the following study 
their phenomenology in full detail, in the SM and beyond. 

In the SM only three basic topologies are present in these decays, depicted in fig. [2J 

• EW penguins 

• CKM- and colour-suppressed tree diagrams 

• Singlet-annihilation diagrams. 

Since the flavour-structure of ir° and p° excludes their production via gluon-exchange, 
annihilation can only contribute if the <j) meson (the flavour singlet) in the final state is produced 
from gluons and the second meson comes from weak (as depicted in fig. [2|) or electromagnetic 
interactions. Since the <f> is colour-neutral and it is odd under charge-conjugation, at least three 
gluons are needed, so that the singlet-annihilation amplitude is formally of higher-order in a s 
and does not appear in QCDF at the next-to- leading order 0{a\) 10]. However, annihilation 



topologies in general do not factorise and cannot be calculated perturbatively, because the 
exchanged gluons may be soft. This means that we can, from a theoretical point of view, only 
rely on the suppression of these contributions by A/mg, where A is a non-perturbative scale, 
and by 1/N C . This leads to the expectation that both the tree and the EW penguin amplitudes 
can receive corrections of 10% — 20% from singlet-annihilation. However, we can also argue 
from a phenomenological point of view that (^-production from three gluons is suppressed by 
the OZI rule 5114541] and should thus be only a small effect, even though this rule is theoretically 
not well understood. In short, our reasoning leads us to the conclusion that in order to test 



NP in B s -> 
uncertainty. 



bit , we have to look for new effects which are much larger than this intrinsic 



In all our 
refs. HEE3- 



calculations of B s — > 4>p , (jm we use the full QCDF decay amplitudes, see 
However, since these are quite involved, we now quote simple approximative 
formulas which can be used as building blocks for an easy calculation of various observables 
such as branching fractions, CP asymmetries and polarisation fractions. Neglecting singlet- 
annihilation we can parameterise the amplitudes in analogy to eq. (|10p as 



V2A(B S -> 4>M 2 



pM 2 1 1 



(19) 
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with Mo, representing a 7r°, a longitudinal p° or a p° with negative helicity. The positive helicity 
amplitude can be neglected in the SM because of its h^ CTl /m 2 B - suppression. We have factored 

out the EW penguin amplitude P EW anticipating its dominance over the colour-suppressed tree 
represented by the tree-to-penguin ratio Tq 2 . A new contribution to the B — > ttK amplitudes 
of the form (fTHj) would also enter the B s — > (j)ir°,(J)p amplitude (fTU|) modifying it as 

V2A(B S -> cj>M 2 ) = P^j (l - e*> + e~ iS ) (20) 



where f E ^ contains a strong phase and 5 is the weak phase introduced in (|13p . If we assume 
the new contribution to be of the order of the SM EW penguin, as required by a solution of the 
"A^cp-puzzle" , we have |f™y| ~ 0(1) and expect a large enhancement of the B s — > <f>TT ,(j)p 
branching fractions, up to an order of magnitude. In order to obtain the same effect within the 
SM one would have to assume an even larger enhancement of the soft non-perturbative physics 
entering the colour-suppressed tree topology in r^ 2 . 

Choosing a phase convention such that -P^v * s rea l; we fi n d 

Pew = 6-45^1 • 10" 9 , P£» = 9-95±§;?l • 10" 9 , = 4.2^ • 10" 9 , (21) 

for the isotriplet meson being -/r , longitudinal p° and p° with negative helicity, respectively. 
We further have 

r TT _ n 41 +0.37 n iq+0.30 ■ 

r c - u.4l_ 041 - u.ld„ 030 i, 

r P,0 _ OO+0.35 _ q ,0+0.28 ■ 
'C — u - oy -0.39 u - lo -0.29 4 ' 

rg" = 0.21^ 9 + 0.15;^z. (22) 



Inserting these numbers into eq. (|19j) we obtain a good approximation of the SM amplitudes 
for the B s decays. Replacing 7 — > —7 in eq. (fT9j) yields the corresponding CP-conjugated 
amplitudes (B s decays). Subsequently one can use the formulas in Appendix IA.4I to convert 
the amplitudes into physical observables. In section [3TT1 we extend these prescriptions to physics 
beyond the SM. 

One should keep in mind that the numbers above are calculated using state-of-the-art 
values for the non-perturbative input parameters, summarised in Appendix IA.31 They are 
based on lattice QCD, QCD sum rules and experimental data. Since our knowledge on these 
parameters is hopefully going to improve in the future it is desirable to have an additional 
parameterisation of the decay amplitudes where the non-perturbative input can be changed. 
We find the dominant sources of theory uncertainties to be (ordered by importance) 

• the form factors ^""^(O) and F± s ~^(0), 

• the CKM angle 7, 

• the non-factorisable spectator-scattering amplitudes, parameterised by the complex num- 
ber Xh and the first inverse moment \b s of the S s -meson light-cone distribution ampli- 
tude. 
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The remaining uncertainties, stemming from decay constants, Gegenbauer moments, quark 
masses and CKM parameters, are much less important so we do not need to display them 
explicitly. Setting the less important theory parameters to their default values we arrive at the 
following approximate expressions for the quantities in eqs. (|21l2'2p : 

P£ W = 17.0^(0) • 10-, rl = -0.12* - 0.02 + ^gg^^ 
B.^*, n \ ln -B p,0 „ 10 - nno , 0.01GeV(l + X H ) 



= 26.2^^(0) • 10~ 9 , rg° = -0.13? - 0.02 + 



p£- = 6.6F B ^(0) • 10" 9 , rff" = 0.14* - 0.06 - O- 02 ^^ 1 gg) . ( 23 ) 
ew - v ; c F B ^(0)A Bs 

The tree topologies r^f 2 suffer from the large spectator-scattering uncertainties due to a strong 
cancellation between the leading order and QCD vertex corrections. Again one can insert these 
formulas into eq. ()19p , this time with arbitrary values and uncertainties for the form factors and 
spectator-scattering parameters, and use the definitions in Appendix I A. 41 to calculate physical 
observables. CP conjugation again amounts to replacement 7 — > —7. 

We conclude this section quoting our QCDF results for the SM values of the B s — > 4>tt°, 4>p° 
observables. As for the CP-averaged branching fractions we obtain 

Bi(B s -> <M°) = l.etH ■ 10" 7 , Br(B s -> <f>p°) = 4.4±g;? • 10~ 7 . (24) 



For comparison we also quote the approximate result according to (|23p : 

Bi(B s -> <jm Q ) = 1.6±H • 10" 7 , Br(B s -> #°) = 4.4±g;| • 10" 7 . (25) 

The smallness of the SM branching ratios compared to other hadronic B decays is due to the 
absence of QCD penguins and non-suppressed tree-level contributions. The measurement of 
these branching fractions is thus challenging and has not been achieved yet. However, we will 
show in later chapters that NP in EW penguins has the chance to enhance the BRs by up to 
an order of magnitude, such that this measurement is a very interesting project. We expect 
that LHCb will be able to measure Br(B s — > 4>p°) while the B s — > <fnr° mode is more suitable 
for a super B-factory where a full reconstruction can cure the notorious difficulties with the 
identification of neutral pions. In case of a strong enhancement B s — > (pp° should also be visible 
in the Tevatron data [55J. The branching ratio Br(i? s — > 4>p°) is dominated by the longitudinal 
polarisation state as can be seen in 

Br(B s -> <f> L p° L ) = 3.7^-5 . l(T 7 (26) 

and the longitudinal polarisation fraction 

h = 0.841°;?* . (27) 

As stated above, one of the main sources of uncertainty in the QCDF predictions is the 
form factor Aq"^. It can in principle be eliminated by considering the ratios 

Bv(B s -> #°) _ . 35 Bt(B s -> 4> L p\ ) _ 010 

By(b s -)■ <m°) " a23 ' By(b s -)■ <m°) " aos ' { ' 
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NP could still be visible in these ratios because in many scenarios it enters B s — > (j)p° and 
B s — > <pTT° in different ways. The cancellation of A^ 3 ^"^ also occurs in the ratios 

Br(-B, ->• <frfl - n nnv+0-008 

Br(B s -> #) " U - UU7 -0-004, 

Br{B s -» _ +Q.023 g£Cgg ->■ <foP° ) _ n m 7 +o.oi9 , 9q ^ 

Br(B a -+ #) " U - U2U -o.oio> Br( J B s -> #) " U - Ui7 ^0.009- l^J 

There however this gain is compensated by additional uncertainties arising from the QCD- 
penguin-dominated decay B s — > (fi<f>. The experimental benefit in these last ratios is that at 
LHCb absolute branching ratios cannot be measured because the absolute number of B s mesons 
is unknown. Finally, we find the direct CP asymmetries to be very uncertain: 

A^(B S -> 4>P°) = 0.19i°;j® A^(B S -> <fm Q ) = 0.27+^. (30) 

Due to the smallness of the branching ratios, these CP asymmetries are also difficult to access 
experimentally, therefore we will not consider them any further. 



3 Model-independent analysis 

In the previous chapter we proposed to test the hypothesis of NP in the EW penguin sector, 
as suggested by the discrepancy in the B — >• ttK observable A^lcp, by a measurement of the 
decays B s — > 07r°,^>/) . In this chapter we support our proposal by a quantitative analysis 
pursuing the following strategy: We parameterise NP in EW penguins in a model-independent 
way by adding corresponding terms to the Wilson coefficients ci'\ ...,C^). By performing a 
X 2 -fit we determine the NP parameters in such a way that they describe well the B — > ttK 
data. In particular they should allow for a solution of the A^4cp discrepancy. Further hadronic 
decays like B — > pK, irK* , pK * are used to impose additional constraints at the 2 a level. 
With respect to the resulting fit we study the decays B s — > (pTT°,(j)p and quantify a potential 
enhancement of their branching fractions. Note that such an exhaustive analysis, correlating 
different hadronic decay modes with sensitivity to isospin violation, is only possible if hadronic 
matrix elements are calculated from first principles like in the framework of QCDF. A method 
based on flavour symmetries, as it has been used in most studies of B — > irK decays so far, 
could not achieve this. In particular, the decays B s — > (pTT°,(j)p , which are our main interest, 
are not related to any other decay via SU(3)f so their branching fractions cannot be predicted 
in this way. 



3.1 Modified EW penguin coefficients 

In the SM the Wilson coefficients CV,...,Cio obey the hierarchy |Cg| 3> |CV| 3> | Cs I > I Cio I a * 
the electroweak scale. This is because Cg receives 1 / sin 2 #yK- en hanced contributions from Z- 
penguin and box diagrams in contrast to CV, while Cs,io are generated for the first time at 
two-loop level due to their colour structure. For our model-independent analysis we consider 
arbitrary NP contributions to the coefficients Cj and Cg as well as to their mirror counterparts 
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Cf p (m 6 )/a 


C? p, (m b )/a 




-0.966 q 7 + 0.009 g 9 


-0.966 g 7 + 0.009 q' 9 


c 8 


-0.387g 7 + 0.002 g 9 


-0.387 g 7 + 0.002 g 9 


Cg 


0.010 g 7 - 1.167ft 


0.010<? 7 - 1.167 q' 9 


Cio 


-0.001 gr + 0.268 qg 


-0.001 q' 7 + 0.268 g 9 



Tabic 2: NLO short-distance coefficients of the EW penguin operators at the scale rrih- Modifications 
to other short-distance coefficients arc negligible. 



C 7 and C 9 . Normalizing the new coefficients to the SM value C 9 ° defined in eq. (|59p in the 
appendix, we have 



cif P (M w ) = C™(M W ) q\% q% = \q%\e^K (31) 



where </> 79 are new weak phases. The coefficient C 9 ° contains the parts of C 9 M enhanced by 
nit/Myy and 1 / sin 2 9y/, as explained in Appendix lA.il There we also describe the scheme which 
we use for the renormalisation-group evolution. Applying it to the NP coefficients leads to the 
low-scale values displayed in tab. [2j They can be compared to the dominant SM coefficient 
Cl M (m b )/a = -1.203. 

In our analysis we will study several different scenarios. First, we consider the cases where 
only one of the coefficients qj, qg, q' 7 , q' 9 is different from zero. This means we assume the 



dominance of an individual NP operator as it has also been done for example in ref. [41]. 
Second, we consider the possibilities of having q 7 = qg, q' 7 = q' 9 , q 7 = q' 9 and q 7 = qg. Finally, 
we study parity-symmetric new contributions corresponding to the three cases qi = q' 7 , qg = q 9 
and q-j = q' 7 = qg = q' 9 . Each of these scenarios can be described by means of two real 
parameters, the absolute value \q\ and phase (j> of the NP contribution under consideration. 
This reduced number of free parameters allows us to perform a fit to B — > nK data and to 
draw meaningful conclusions on the B s — > (j)7r°,(j)p decays. The study of this large set of 
well-motivated simplified scenarios is assumed to represent all relevant features of the general 
framework with unrelated q?, qg, q 7 , q 9 . 

Our main motivation for adding NP to the coefficients C 7 \C 9 ^ was the claim that the 
AAcp discrepancy can be solved in this way, namely by generating the terms tewi ^ew 
fg W introduced in eq. ()13|) . Introducing individual terms for each of the four relevant Wilson 
coefficients, we obtain 

E nw.i*-* = (<* -4) [(-0- 12 )±ao5 + (-0.02)i^i] + 



1=7,9,7' 9' 



(q 9 -q'g) [0.12±g;gS + 0.021°:™*] , 
E = (97- (6) [0.10i^ + 0-Oll^i] + 

(qg-q'g) [0.04l° :o f + (-0.005)1^26 



j=7,9,7',9' 
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£ r 1w,i^ = (QT-4r) [0-03 



+0.04 
-0.07 



+ (-o.o6)±8:S?i] + 



«=7,9,7',9' 



(?9 



g 9 ) [0.007; 



•+0.003 
-0.010 



+ (-0.006) 



+0.012 - 
-0.003 4 



(32) 



Let us briefly discuss the main characteristics of these coefficients: 

• First of all, note that parity-symmetric models obviously do not contribute to B — > ttK 
at all. This general feature of B decays into two pseudoscalar mesons (PP decays) follows 
from eq. (|63|) . Therefore such a scenario cannot solve the A^Lcp discrepancy. 

• The contributions f EW 7 (,) and f EW 9 (/) tend to cancel each other. Hence in the scenarios 
with qj = qg and q' 7 = q ! 9 only a negligible new colour-allowed EW penguin contribution 
is generated. 

• Whereas Re(f^ w g(/) ) features the typical colour-suppression with respect to Re(f EW g (/)), 

this pattern is not obeyed by the q? terms. This is due to a conspirative interplay of 
the large mixing of Cy into (compare tab. [2]), constructive interference of the new 
ci and Cg^ contributions in the QCDF coefficient Og and a chiral enhancement factor 
r^ K « 1.5 multiplying in eq. (|6ip for the topological amplitude. None of these three 
effects is present in the q§ case. 

• The annihilation coefficient {/) develops a large imaginary part. In scenarios with 
non- vanishing qlP this term gives the dominant contribution to AAcp. 

From eq. f|15[) we see that the A^4cp discrepancy can be solved either through tew or 
through r EW . Except for the parity-symmetric models, all the scenarios mentioned above can 
achieve such a solution. In fig. [3] this is illustrated for the cases with a single q? or q$ and 
for the qj = qg scenario. Graphs for the respective mirror scenarios are obtained by a 180° 
rotation. The yellow region contains those points of the (Re(gj), Im(gj)) - plane for which the 
theory error band overlaps with the experimental 1 a region, whereas the blue region represents 
those points for which also the experimental central value lies within the theory error interval. 
The red circle illustrates the minimal \q\ - value needed to reduce the A^4cp tension below the 
lo" level. For the three scenarios in fig.[3]we read off \q-j\ > 0.3, \qg\ > 0.8 and \q-?\ = \qg\ > 0.4. 
The fact that in the qj = qg case only a small NP contribution is needed, in spite of the absence 
of teWj demonstrates the importance of the annihilation term f EW . Finally, we like to stress 
that the solution of the AAqp discrepancy via a minimal \q\- value requires the adjustment 
of the phase to a certain value. Realistic scenarios avoiding such a fine-tuning have larger 
| q\ -values, typically |g| ~ 1. 

Our main goal is to study the impact of such a NP scenario on the decays B s — > (fin , (fip°. 
The NP contributions to , . . . , c{q generate the fj^ - terms introduced in eq. (|20|) . Intro- 
ducing again individual terms for contributions from the various Wilson coefficients, they read 
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Figure 3: NP contribution needed to solve the AAcp discrepancy in the three scenarios (from left to 
right) with single q 7 , single qg and equal qj = qg contribution. Yellow region: Theory error band and 
experimental 1 a region overlap. Blue region: Theory error band and experimental central value overlap. 
Red circle: Minimal magnitude of the NP contribution needed to reduce the AAcp discrepancy below 
the 1 a level. 



for the four different amplitudes 



i=7,9,7',9' 



E 



~P,0 -i8i 
' EW, i e 



i=7,9,7',9' 

E 

i=7,9,7',9' 



' EW, i e 



V fP<+ P ~iSi 
i=7,9,7', 9' 



-0.9 (q 7 + q 7 - q 9 - q' 9 ) , 



0.9 (q 7 - q' 7 + q 9 - q' 9 ) , 



-0.6 (q 7 + qg) , 



0.6 (q' 7 + q'g) X P^W/^EW 



(33) 



where we have neglected c^g-contributions to fg W and q 7 9 -contributions to fg W according 
to their Aq CD /?ti^ suppression. The SM EW penguin amplitude drops out of the total 

expression ([20]) of the amplitude, Pgfc is given in eq. (|21 j) . The parameters f E -^ develop only 
very small strong phases and uncertainties not indicated in (|33[) . This is because they are ratios 
of equal topologies such that uncertainties and strong phases approximately cancel. 



We have stated the expressions in eq. (j33[) for two reasons: Firstly we want to show the main 
consequences of non-vanishing qig for the B s — > (/>7r°,(/>p decays. We see that for \qi\ = 0(1) 
indeed new contributions with the magnitude of the leading SM EW penguin are generated. 



57T 



While parity-symmetric NP was invisible in B — > ttK, it could be detected in B s 
and in principle also in B s — > (f>p° due to the different interference patterns of and 
with the corresponding SM contributions. Furthermore, left- and right-handed NP could be 
distinguished by a polarisation measurement of B s — > (j)p°. This general feature of decays to 



vector- vector final states has been pointed out by Kagan 56J]. Note that the question of left- 
vs. right-handed NP cannot be answered from B —> irK alone since, as we have seen, the two 
scenarios differ only by a rotation in the NP parameter space. 



The second benefit of eq. (|33|) is that it allows for a simple calculation of B s 



?7T, 
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observables to a very good accuracy. In this way it permits a study of these decays without 
the extensive implementation of QCDF. One simply evaluates the amplitude eq. (|2Up inserting 
eqs. (|21|22p - or alternatively eq. (f23j) - and the NP part from eq. ([33]) . The CP-conjugated 
amplitude is obtained by flipping the sign of 7 and replacing the qi by their complex conjugates. 
Subsequently one can use the formulas in Appendix IA. 41 to calculate observables. 

3.2 Fit to B — y nK data and constraints from other decays 

The four B — y irK channels are the most precisely measured hadronic b —y s decays. For this 
reason, we use experimental information from these channels as input for our quantitative NP 
analysis by performing a fit of qtj) and qg to B — y irK data. This fit will be an indication for 
values of the g- that are realistic to expect and consequently will enable us to find an expected 
size of enhancement of the branching ratios of B s —y (f)p° and B s — > 07r°. In the treatment of 
theoretical and experimental uncertainties in the fit we follow the RG.t scheme [?J ■ More details 
on this issue are given in Appendix [Cj 

Since B — y irK decays are greatly dominated by QCD penguins and since they suffer 
from large theoretical uncertainties, it is obvious that NP effects residing in EW penguins are 
difficult to find in branching ratios and direct CP asymmetries. It is more useful to consider 
instead particular combinations of these basic observables which highlight the isospin- violating 
contributions and can be predicted with a better precision. For example, it is difficult to see 
a need for isospin-violating NP by looking at the two CP asymmetries entering in A^4cp ( see 
tab.[T]) since both of them have more than 100% theoretical uncertainty, reflecting the difficulty 
of predicting strong phases in QCDF. In the difference A^cP) the theory uncertainties cancel to 
a large extent such that the discrepancy with experimental data becomes clearer. Moreover, in 
realistic models a new contribution in the EW penguin sector usually comes in combination with 
NP of comparable size in the QCD penguins since the new contribution in general matches onto 
a linear combination of the QCD and EW penguin operators. By considering mainly isospin- 
violating observables, one reduces the sensitivity to additional effects from new contributions 
to C3, ...,Cq which we did not include in our parameterisation (|3ip . 

For our B —y irK fit, we use suitable ratios of branching fractions and differences of CP 
asymmetries as well as the mixing-induced CP asymmetry in B° —y ir°K as input. The 
definitions of these quantities are summarised in Appendix [Bj Many of these quantities have 
also been considered in the past in the context of flavour symmetry analyses of B — y ttK decays. 
A summary of theoretical predictions vs. experimental results for all these observables as well 
as for the B —y irK branching fractions and CP asymmetries is provided in tab. [TJ At present 
none of these quantities deviates from the SM prediction by more than la apart from A^4cp- 
We thus expect AAqp (and to a lower degree also Sqp) to pull the fit towards non-zero values 
of the q^ whereas the other observables will favour values close to the origin of the complex 
plane. 

In addition to the fit we consider constraints on the NP parameters arising from a large 
number of hadronic B decays. To this end we compare the theoretical prediction of an observ- 
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able as a function of the g- to its measured value and extract a 2c-constraint as follows: 



(^-theo ~t~ ^theo, sup) ^ (^cxp 2(7 cxP)U1 f) 

Point [q^' ; 9 space] = <( [ and (x thco - cr thc0iinf ) < (x cxp + 2cr cxPjSup ), (34) 

excluded otherwise. 



(/) allowed if 



Here (xtheo)-^!' 50 'ZJ re P resen t s the theoretical prediction for the respective physical observable. 
The uncertainty does not imply a particular probability distribution but the true value is 
supposed to lie within the error interval. For the experimental value x exp ± cr exp a Gaussian 
error is assumed. 

This procedure is applied to data from B — > ttK as well as to data from the B — > pK, 
B — > ttK* and B — > pK^ decay channels, which are simply the pseudoscalar-vector (PV) and 
vector-vector (VV) modes corresponding to B — > ttK. The PV and VV modes turn out to 
be more sensitive to isospin-violating flavour topologies than their PP counterparts because 
the leading QCD penguin amplitude is smaller. Experimental information on these decays, 
however, is not (yet) as precise as the available data for the irK modes. Therefore we do not 
include the PV, VV modes into the fit but prefer to consider them as constraints at the la 
level only. Nonetheless, the constraints from B — > pK and B — > ttK* give some information 
complementary to the one from B — > ttK because they test different chirality structures than 
B — > irK and are therefore sensitive to other linear combinations of the g^'g. Moreover, we 
apply eq. (|34p also to data from B — > K^4>, B s — > (fxj) and B s — > KK decays even if they only 
carry a small sensitivity to EW penguins. 



3.3 Results of the model-independent analysis 

We now discuss the results of the analysis outlined in the previous section. The aim is to make 
predictions for the B s decays in combination with the regions of the g|' g parameter space which 
are preferred, or not yet excluded, by experimental data from B — > ttK and related decays. 

In figs. H] and [5] we present graphs showing the enhancement Br SM+NP /Br SM of the 
B s —7- (f)p°,(fiir branching ratios as a function of the parameters qp in some representa- 
tive scenarios. The displayed numbers are obtained with our default hadronic input. In order 
to be distinguishable from the SM, a particular scenario must at least provide a value for 
Br SM+NP /Br SM which exceeds a potential enhancement factor faked by hadronic uncertainties 
in the SM prediction. Parameter points for which the enhancement factor lies within the the- 
ory error band of the SM prediction are represented by the red-striped ring. The SM itself 
corresponds, of course, to the origin of the plots and is highlighted by a black dot. 

In our sample models we introduced NP exclusively in the EW penguin operators. In 
realistic models, however, a new contribution in the EW penguin sector usually comes in 
combination with NP of comparable size in the QCD penguins since the new contribution in 
general matches onto a linear combination of the QCD and EW penguin operators. In order 
to account for this fact, we use the experimental data in two different ways: 

a) We present a fit using isospin-sensitive quantities in the B — > ttK decays, such as ratios 
of branching fractions and the differences of CP asymmetries as discussed in Appendix [B] 
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bcenano 


Br^-x/m ) 


Br{B s -+<j> L p u L ) 


Br(B s ^<t>p°) 


Br SM (B s -^7r°) 


Br aM (B S ^4>LP° L ) 


Br BM (_B s -></>p°) 


q 7 ^0 


11.0 (18.7) 


6.0 (9.9) 


5.3 (8.4) 


99^0 


8.8 (16.4) 


9.3 (17.0) 


8.7 (15.1) 


q 7 = q 9 / 


1.0 (1.7) 


11.5 (21.1) 


10.8 (18.7) 




8.3 (15.6) 


8.8 (16.4) 


8.4 (14.7) 


9 9 /0 


6.2 (9.8) 


2.8 (5.6) 


2.7 (5.0) 


57 = 9 9 ^ 


1.0 (1.7) 


5.7 (8.9) 


5.2 (7.9) 


97 = 99 + 


1.0 (1.7) 


1.0 (1.7) 


1.6 (2.3) 


g 7 = g 9 7^ 


1.0 (1.7) 


1.0 (1.7) 


1.1 (1.8) 


97 = 9 7 + 


29.5 (48.1) 


1.0 (1.7) 


2.1 (3.0) 


99 = 9g / 


11.1 (20.5) 


1.0 (1.7) 


1.5 (2.2) 


97 = 97 = 99 = 9 9 + 


1.0 (1.8) 


1.0 (1.7) 


2.3 (3.4) 



Tabic 3: Maximal possible enhancement of branching ratios compatible with the constraints from B — > 
wK, pK,irK* , pK* ,4>K,4>K* and B s — > cjxf),KK decays at the 2a -level as well as the ler confidence 
level from B —> irK decays (method a)). We use the default SM value and the default (the maximal 
theoretical) total value of each branching ratio for an optimally chosen qi value. Numbers in the second 
part of the table are obtained ignoring A^4cp • 

and £3 plus the time dependent CP asymmetry Scp{B° — > irK ). In figs. I4|5I we indi- 
viduate the la region by a solid black line. At the same time we consider the constraints 
one obtains from all non-leptonic B — > ttK, pK,irK* , pK* , <j)K, <j)K* and B s — > cf)(f>,KK 
decays at the la - level and mark the allowed region by a (dark-)green area. 

b) We exclude the time dependent CP asymmetry Scp(B° — > ttK ) from the fit and we 
restrict the constraints to the subset of observables which are particularly sensitive to 
isospin violation, see Appendix|Bj This procedure enlarges the la confidence-level of the 
B — > ttK fit (indicated by the grey line) as well as the region of 2c-allowed parameter 
points by the areas depicted by lighter colours in the figures. 

Whereas the results from a) are valid only if NP is strictly limited to the electroweak penguin 
operators, the results from b) are expected to remain approximately valid also in presence of 
NP in QCD penguins, since such a kind of NP has only a minor impact on the quantities 
considered in b). 

We find that the B — > rrK and related decays set quite strong constraints on the parameter 
space, especially in scenarios where qg 7^ or q' 9 / 0. This basically rules out the possibility 
of having > 5, i.e. NP corrections cannot be much larger than the EW penguins of the 
SM. The fact that the SM point is always excluded at the 2 a level is a direct consequence of 
the AAcp data. According to the sign pattern in eq. (|3"2"j) . the B — > ttK fits of the primed 
and unprimed scenarios in figs. 14151 are related to each other through rotation by 180°. The fit 
works best in the t/g^ scenario where (using method a)) the best fit point is given by 

\$\ = 1.9 <$ = -100° (+180°). (35) 
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bcenario 


Br^-x/m ) 


Br(_B s — s- 




Br(B s 




Br SM (B s -^7r°) 


Br SM (5 s - 




Br SM (B s ^4,p°) 


q 7 ^0 


77.4 (134.1) 


72.5 ( 


117.6) 


66.9 ( 


104.7) 


99^0 


12.0 (21.9) 


12.6 


(22.8) 


11.8 


(20.3) 


97 = 99 + 


1.0 (1.7) 


52.9 


(90.9) 


49.4 


(81.0) 




56.6 (99.2) 


59.5 ( 


103.2) 


54.0 


(90.5) 




13.0 (20.6) 


13.0 


(20.5) 


11.7 


(18.1) 


57 = 99 ^ 


1.0 (1.8) 


36.3 


(58.2) 


32.8 


(51.2) 


97 = 9 9 / 


1.0 (1.7) 


1.0 


(1.7) 


2.5 


(3.8) 


q' 7 = q g ^0 


1.0 (0.0) 


1.1 


(1.8) 


2.5 


(3.8) 


q 7 = q 7 ^0 


76.0 (131.9) 


1.0 


(1.7) 


3.8 


(5.5) 


99 = q' 9 / o 


13.0 (20.6) 


1.0 


(1.7) 


1.5 


(2.2) 


97 = <? 7 = 99 = 9 9 + 


1.0 (1.8) 


1.0 


(1.7) 


4.0 


(5.9) 



Tabic 4: Maximal possible enhancement of branching ratios compatible with the constraints from isospin 
sensitive observables in B — >• irK, pK, irK* decays at the 2a - level as well as the 1-er confidence level 
from B — > wK decays, without including Scp(B° — > ttK°) (method b)). We use the default SM value 
and the default (the maximal theoretical) total value of each branching ratio for an optimally chosen qt 
value. Numbers in the second part of the table are obtained ignoring AAqp- 



This parameter point yields a full agreement of all the B — > ttK observables with the experi- 
mental mean values (for Sqp the nearly exact value £*cp — 0.55 is obtained). In the q 7 ^ = q^ 
case a plateau of x 2 = points arises due to the large theoretical errors. It turns out that 
the B —7- ttK observables are not very sensitive to the g^-only scenarios and so the fit does 
not work well here. Hence within the g^-only setting one can mainly rely on the 2 a con- 
straints. The total sets of constraints of the primed and unprimed scenarios are not anymore 
related to each other in a simple way since they involve PP decays subject to a 180° rotation 
together with PV decays which are unaffected by a q 7 $ «->■ q' 79 exchange. It turns out that 
the constraints are stronger in the (/g-only and in the q 7 = q' 9 scenarios than in their unprimed 
counterparts and that the best fit regions are cut away in these cases. 

Prom figs. I4|5l the enhancement Br SM+NP /Br SM of the B s branching fractions can be read 
off with respect to the different constraint- and fit-regions. A large enhancement of the 
B s — > <j)p°,(j)7{ branching ratio can be expected in many scenarios, especially in those in- 
volving q 7 ^ 0. The fact that large parts of the allowed regions do not overlap with the SM 
uncertainty regions is encouraging. It means that, if such NP is realised in nature, it could 
be possible to probe it easily. In tabs. [3] and H] we quote the maximal enhancement factors 
that can be obtained considering all points in parameter space which lie within the 1 a region 
of the B — > ttK fit and fulfill the additional 2 a constraints. The numbers in tab. [3] refer to 
input a) while the numbers in tab. prefer to input b). The first number in each cell represents 
Br^J l " NP /Br^ d while the number in brackets represents Br| 1 I ^ hNP /Br| 1 ^ [ d , both evaluated for 
the qi value which gives the largest enhancement. Here "max" and "med" refer to the upper 
limit of the theoretical uncertainty range and to our default value, respectively, according to 
our input given in Appendix IA.31 Exploiting the theory error in favour of an enhancement, 
the number in brackets gives the absolutely maximal enhancement possible for each scenario 
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Observable 


M = 1.9, 
09 = -100°, 


Br(B s -> (fin ) ■ 10 6 
Br L (B s -> • 10 6 
Br(B s -> #°) • 10 6 


u - oo -0.19 

n qn+ - 98 

1 1O + 0.95 

1 - LO -0.38 



Table 5: Values of various obscrvablcs at our best fit point in the scenario qg ^ 0. 

whereas the first number gives a typical enhancement factor, but still for the most enhancing 
parameter point. Concerning the parity-symmetric scenarios one should have in mind that like 
the SM they violate A^4cp at the > 2 a level since they have no impact on B — > nK decays. 
The corresponding enhancement factors shown in tabs. [3] and H] are obtained ignoring A^4cp 
but taking into account all other constraints. 

In most scenarios an enhancement of more than an order of magnitude is possible. Ex- 
ceptions are B s — > (fin for q 7 = q§ and B s — > 4>lPl f° r parity-symmetric models and have 
their origin in the pattern of eq. (|33p . Furthermore, effects in the q'g and the q' 7 = q'g scenarios 
are limited by the small allowed region resulting from the B — > nK fit. Largest effects oc- 
cur as expected in the scenarios which are least constrained by B — > nK, i.e. the single q^f 1 
and the parity-symmetric models. Especially in these cases a B s — > (fin measurement would 
complement B — > ttK data and, while the parity-symmetric models lack the motivation via 
the A^Lcp discrepancy, the q' 7 setting resolves it with ease (see Fig. [3]). Moreover, we like to 
stress that B — > irK data alone cannot distinguish among opposite-parity scenarios because 
such scenarios generate equal results for the B — > ttK observables (for 180°-rotated parameter 
points). Therefore an analysis of B — > ttK should for example be supported by the analysis of 
a PV decay, suggesting B s — > (fin as an ideal candidate. 

We have seen that NP in the EW penguin coefficients allows for an enhancement of Br(£? s — > 
(fin, <fip) of more than an order of magnitude. According to the simple topological structure of 
these decays, the observation of such an effect would be a clear and unambiguous signal for such 
a scenario. It is interesting to raise also the reversed question, i.e. whether the absence of such 
an effect would rule out a NP solution of the AAqp discrepancy, at least for a specific scenario. 
This is, however, not compulsory. In nearly all the considered settings there are points within 
the 1 a region of the B — > irK fit which do not generate an enhancement of Br(i? s — > (fin, <fip). 
The only exception is the gg-only case: Here an enhancement factor of at least 2.1 would occur 
in B s — > (fin. This time we have exploited the theoretical error in disfavour of an enhancement 
(for the default value the factor is 2.7). Finally we provide in tab. [5]the B s — > (fin, <fip branching 
ratios for the best fit point in the gg-only scenario. 
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Figure 4: Enhancement factors of the B s — > <fip°, 4>ir° branching ratios with respect to their SM values. 
The black dot represents the SM result while the red striped region shows the theoretical uncertainty in 
the SM. The dark green area is the region allowed by the 2 a constraints from B — > tyK^*\ pK^*\ 4>K^ 
and B s (jxf>, KK decays; for comparison, the light green area represents the area allowed by constraints 
from isospin-sensitive observables only, considering only B — > ttK,ttK^*\ pK decays. The solid black 
line represents the la CL of the fit with Scp(B° ttK°), while the solid grey line represents the la 
CL of the fit without it. Here the scenarios (77 7^ (upper row) and qg ^ (lower row) are displayed. 
For the (77-scenario the la region of the fit is the region to the left of the black (grey) curve. 
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Figure 5: Enhancement factors of the B s — > <fip°i 7r ° branching ratios with respect to their SM values. 
The meaning of the contours and regions is the same as in fig. 01 Here the scenarios q' 7 ^ (upper row) 
and <7g ^ (center row) and qi — q$ and q' 7 = q' 9 ^ (lower row) are displayed. For the g^-scenario 
the la region of the fit is the region to the right of the black (grey) curve. 
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4 Analysis of viable New-Physics models 



In view of the results in chapter [31 the question arises which concrete models for NP can provide 
a large new EW penguin amplitude without being excluded by present data. In this section 
we consider a number of well-motivated NP models. The main difference with respect to the 
model-independent analysis is the possibility of adding constraints from other flavour processes 
beyond the hadronic B modes, e.g. the semileptonic decay B — > X s e + e~, the radiative decay 
B —7- X s j and B s -B s mixing. These processes usually yield tight constraints on new flavour 
structures and it has to be investigated if the effects in B — > nK and B s — > (pp°,cpTT survive 
these constraints. 



4.1 Constraints from semileptonic decays and B s — B s mixing 

Before turning to the models we summarise here how we implement constraints from semilep- 
tonic B decays and B s -B s mixing. The inclusive semileptonic decay B — > X s e + e~ is generated 
by electroweak interactions and therefore correlated to the hadronic EW penguins in many 
models. We describe it by the effective Hamiltonian ©, adding the operators 

Qw = (s a b a )v-A (ll)v and Q 10 a = (s a b a ) V -A (11) A (36) 
and the corresponding mirror cop ies Q'gy , Q'inji- The SM expressions for the Wilson coefficients 



can be found e.g. in refs. 57j, [58|]. Following [571 ] and extending the formulae therein to include 



effects of the mirror operators, we use the effective Hamiltonian to calculate the ratio 

■Avib^ se+e-) 
Wtf) ^' r( ^ cep) ■ (37) 

where q 2 = (p e + + p e -) 2 is the squared invariant mass of the lepton pair. This ratio has the 
advantage that its theoretical uncertainty is considerably reduced with respect to the simple 
branching fraction. We integrate over a continuum region below the tp resonances to find the 
integrated ratio 

6GcV 2 

R e+e -(q 2 )dq 2 , (38) 

lGeV 2 



-n-e+e- I [1,6] 

which we can finally compare to the experimental result 



Br e+e -| [1:6] = (1.60 ±0.51)- 10~ b , (39) 

also normalized to the semileptonic decay. We require -R e +e-|[i,6] to be compatible with exper- 
imental data according to (pT 



Besides the inclusive B — > X^e + e~ , also the exclusive mode B — > K*l + l~ has been found to 
be a useful constraint for NP 63H65I]. Here we focus only on the forward-backward asymmetry 
A-p-Q of this process which gives a constraint complementary to that of R e + e - |[i,6l- I n the 
light of present experimental data we require the sign of Afb(q 2 ) integrated over q 2 > 14 GeV 2 
to be negative. 

For completeness we note that we use a renormalisation-group evolution analogous to the 
one of the EW penguin operators, treating the parts of Cgy and C\qa enhanced by xtw = 
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rrit/Myy and/or 1/sin 2 9y/ as leading order. This results in the following SM initial conditions 
at the scale fi ~ 0(M\y): 



a 



(0) 
9 V 



Y (x tw ) _ x t w 
sin 2 9w 2 



a 



in 

9V 



a 

^ f-4Z (x w ) + — + - 



a a s (Yi(x t w) 



C 



(o) 

1CL4 



2tt sin 2 



27r 47r \ sin 2 #^ 
a a s Yi(xtw) 



^x t w [ i ~y 



2vr 4vr sin 2 6> 



ir 



(40) 



The functions Yq \ and Zq can be found e.g. in 661 ]. 



Finally, we consider constraints coming from B s -B s mixing, which is described by the 
effective weak Hamiltonian 



H 



(2) 
off 



4tt 2 



(41) 



with the operators [67 1 

Q VLL = (a a -fP L b a ) (WvPLbf)), 

Ql LL = (s a P L b a ) (spP L bp), 
QLR = ^ a -fP L b a ) (sp^PRbp), 



n SLL 
V2 



Q 



LK 



(s a a^ u P L b a ) {spa^ u P L bp), 
(s a P L b a ) (spP R bp) 



(42) 



and the mirror copies Q VRR , QSRR and Q| R R - in the SM only C VLL ^ 0, while in extensions of 
the SM all operators can receive contributions. The matrix element relevant for B s -B s mixing, 



M 



12 



2m B: 



■(B° S \H%\B° S ) 



is evaluated using lattice results from ref. 



AM, = 21 M 



Besides the B s -B s mass difference 
g.| ^ (17.77 ± 0.12) ps" 1 , 



69J, we use the quantity [7_ 



M 



B S ,SM 
12 



M 



lAJe 



i4s 



(43) 



(44) 



(45) 



as additional constraint. This observable has been analysed in ref. 7l|] in different generic NP 
scenarios and evidence for a NP contribution with a large new weak phase has been found. A 
fit of A s and the analogous quantity A^ to data shows a 3.6 a discrepancy for the SM value 
A s = 1. In our study of the Z' models we take those points of the NP parameter space as 
excluded which give a A s outside the 2 a region drawn in fig. 9 of ref. 7l[ . 
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4.2 The modified- Z°-penguin scenario 

The simplest class of models with large new contributions to EW penguins comprises models 
with a modified Zsb coupling. Such a FCNC coupling can either be generated by integrating 
out new heavy particles, e.g. in super symmetric models or fourth-generation models, or it can 
exist at tree-level in more exotic scenarios like models with non-sequential quarks, see e.g. 



ref. 72J. Consequences for hadronic B decays have been considered for example in [731 ] . more 
detailed analyses of the motivation and the effects in flavour physics have been performed in 
10. 



4.2.1 Effective Theory 



Our parameterisation of the Zsb coupling follows ref. 73j. At the electroweak scale we have 
an effective theory with the Lagrangian 

^ = "Zdb- ^ * - 75) + 4V(1 + 75)] djZ„ (46) 

4 cos t>w 

where /, J are generation indices. Since the flavour-violating couplings are expected to be 
small, the flavour-diagonal couplings of the Z bosons are to a first approximation the same 
as in the SM. Matching tree-level diagrams with Z exchange onto the AB = AS = 1 effec- 
tive Hamiltonian adds new contributions 5C{ to the SM Wilson coefficients Cj and generates 
coefficients C[ of the mirror operators. The resulting contributions at the electroweak scale 
read 

xn — L n* — K R 

sc 7 = \^-We w , c> 7 = coJ e w , 

3 A W 3 X f> 

dc 9 = - 7 \^co S 2 e w , c' 9 = \^Lwe w . (47) 

3 xf> 3 x\ s) 

They reach the size of the dominant SM Wilson coefficient Cg(fiw) if 



.sb 



a 



\Klr\ ~ l«SM| = . 2a \ { t s) C (x tw ) ~ 0.00035, (48) 
7r sin t>w 



where Cq(x) is a loop function, see e.g. ref. [66J. Such a scenario corresponds to ~ 0(1) 
in our model-independent analysis, thus we expect significant effects in hadronic B decays for 
such values of k s £ r . 

From the Lagrangian (|46|) and the SM coupling of the Z to leptons we also obtain corrections 
to the short-distance coefficients of the semileptonic operators (|36p . namely 

5C 9V = - ^ (2sm 2 ^ - IV C ' w = (.2 S m 2 e w - \ 



.sb 



1\ „, itiS \ 
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This enables us to study constraints on ftf 6 R from semileptonic B decays as indicated in the 
previous section. 

Diagrams with Z-exchange contribute also to B s -B s mixing via the Wilson coefficients 
4vr 2 / K s r b \ 2 4tt 2 



V2G F M^ \\<f) J ' V2G F M* 




87T 2 K s T b K s t b 



C 1 LR = °" , ^fWv ( 5 °) 

Explaining the discrepancy in A s defined in eq. ()45p with the help of these new contributions 
would push the couplings k s £ r to large values. Note, however, that in most realistic cases the 
couplings k s l R are loop-induced with the consequence of eq. (|50p actually representing two-loop 
effects. Usually such scenarios provide also one-loop contributions from box diagrams which 
then are more likely to account for the A s discrepancy. Therefore we prefer not to include A s 
as a constraint into our analysis and regard a potential relaxation of the A s discrepancy only 
as a bonus feature. 



4.2.2 Results 

In our study of the modified Z coupling we consider the three special cases of non-vanishing 
K sb only, only and K sb = similarly to the model- independent analysis. Since cos 2 6\y 3> 
s'm 2 9w, the k s £ scenario shares its most important features with the qg setup of the model- 
independent study and the same holds for k r and q' 7 . This expectation is confirmed by the 
graphs in fig. we only note that we get a 180° rotation due to the signs of 5Cg and C' 7 . We 
have again marked the la - region of the B — > irK fit by a black line (as well as the additional 
3a- black dotted line) and the region allowed by the 2 a constraints from all hadronic decay 
observables by a green area. The displayed regions refer to input a) , as defined in section 13. 3| 
while we have refrained from showing the corresponding regions for input b). 

The main difference to the more general model-independent approach is that we now face 
additional constraints from semileptonic decays and B s -B s mixing. The allowed region for the 
former is given by the interior of the blue dashed curve, the allowed region for the latter by 
the orange areas outside the zone preferred by the B — > ttK fit. We see that the A s anomaly 
of B s -B s mixing cannot be resolved in a modified Z scenario when fulfilling at the same time 



the semileptonic constraints. This has already been noted in Ref. [75J. Here we recognise 
that also B — > irK, pK, irK* data are not compatible with a solution of A s in this way. In 
the previous section we remarked that it is plausible to assign the explanation of A s to other 
effects not directly related to the modified Z coupling. Pursuing this strategy, we are left with 
the semileptonic decays which are compatible with the 1 a region of the B — > irK fit for all 
three cases but constrain the FCNC couplings K sb R to very small values as can be seen from 
fig. [6] where the coupling k s £ r is normalised to |ksm| in Eq- ([1BD - 

As a consequence we expect no significant effects in B s — > 07r° , <pp° . This expectation 
is confirmed by fig. [6] and by the maximum enhancement factors given in tab. [6l which are 
determined in analogy to the ones in tab. [3l In the n sb = k r case no enhancement occurs at all 
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Scenario 


Br(B a ^cfm u ) 


Br L (B 3 ^0p u ) 


Br(B a ^</>p u ) 


Br SM (B s ^07r°) 


Br£ M (B a ->&>°) 


Br MM {B a ^<j>p°) 


„sb _ ^sb 
K L — h R 


10.3 (19.0) 
48.3 (84.6) 
1.0 (1.7) 


3.6 (7.0) 
15.5 (28.2) 
1.0 (1.7) 


3.4 (6.3) 
14.2 (24.8) 
1.2 (1.8) 


With additional semileptonic B decay constraints 


„sb _ ^sb 
K L — K R 


1.6 (3.0) 
4.0 (6.5) 
1.0 (1.7) 


1.1 (2.2) 
2.4 (3.9) 
1.0 (1.7) 


1.1 (2.0) 

2.2 (3.5) 
1.1 (1.7) 



Tabic 6: Maximal possible enhancement of the B s branching ratios in the modihed-Z°-pcnguin model. 
The upper part has been calculated with the method a) of the model-independent analysis, the lower 
part includes the 2<r constraints from semileptonic decays. 

because of the pattern in eq. (j33|) : Equal contributions to C7 and Cg and to Cg and Cj cancel 
pairwise. The largest effect which one could gain in the other scenarios is a factor of ~ 4 in the 
case where only 7^ 0. Therefore an enhancement of B s — > 07r°, (ftp due to a new modified 
Z contribution in practice becomes indistinguishable from the potential enhancement caused 
by a large non-factorisable SM effect. In fig. [6] this is reflected by the fact that the red-striped 
ring representing parameter points which reproduce the SM result for the B s decays nearly fills 
the whole allowed region of the parameter space. 

Our results can be summarised as follows: The constraints from semileptonic decays still 
allow for a solution of A^4cp via a modified Z coupling. This possibility would be excluded if 
an enhancement of B s — > cjm° or B s — > (ftp by an order of magnitude was found. 
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Figure 6: Enhancement factor for the B s — > tfip ,(j)ir branching ratios with respect to their SM 
values in the modificd-Z°-pcnguin scenario. The green area represents the region allowed by the 2a 
constraints from all the considered hadronic decays, while the area inside the dashed blue line represents 
the region allowed by the 2er constraint from semi-leptonic decays. The areas inside the dashed orange 
line represent the parameter values for which the modificd-Z -penguin would solve A s . See the text for 
further explanations. 
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4.3 Models with an additional U(l)' gauge symmetry 



The presence of a heavy Z' boson associated with an additional U(l)' gauge symmetry is a 
well- motivated extension of the SM. This additional symmetry has not been invented to solve a 
particular problem of the SM, but rather occurs as a byproduct in many models like e.g. Grand 
Unified Theories, various models of dynamical symmetry breaking and Little-Higgs models. In 
many scenarios the Z' boson is expected to have a mass at the TeV scale. It also appears in 
the form of a Kaluza-Klein excitation of the SM Z° boson in theories with extra dimensions. 
An extensive review about the physics of Z' gauge-bosons can be found in fl$\ . Here we are 
interested in implications for flavour physics as discussed in 73, 77-7^]. 



4.3.1 Effective Theory 

We consider a model with an additional Z' gauge-boson, neglecting Z-Z' mixing and assuming 
the absence of exotic fermions which could mix with the SM fermions through non-universal 
Z' couplings. We write the general quark-antiquark-Z' coupling as 73. l79| 



reff 



9U(1) 



I J 



iV(i 



75 J 



+ CflV(l + 7&)] djZ'r 



(51) 



and similarly for the up-type quarks. The couplings of interest are the flavour-changing Clr 
as well as the flavour-conserving charges (1 R = Q^ U R and Clr = Clb- Note that SU(2)i 
invariance implies Cl = Cl = Cf whereas no restrictions hold in case of C R . Following 
ref. |73| we introduce the parameter 



9u(iy 



Ml 



(52) 



with gjjny denoting the gauge coupling of the additional U(l)' gauge group and Mz' being the 
mass of the Z'-boson. We then find the following additional contributions to the short-distance 
coefficients at the electroweak scale: 



5Ca 



sc 7 
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3 is-) 

0. 



9 /^sfe 



c> 



c> 



3A 



sb 
R_ 
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/-sb 
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R 



(53) 



Apart from ^, FCNC transitions are controlled by the free parameters Cl J r- Depending on them, 
the flavour-changing transitions contribute in general to both QCD and EW penguin operators, 



as well as to their mirror copies. Here we follow the approach of refs. [45j, I77H79I] in which the 



main contribution is supposed to reside in the EW penguins, i.e. \5C^^{nw)\ *C \5Cj(nw) 
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1^3,5 (a*w) I ^ I^yCmw)!- We implement this assumption by setting + 2^ = (1 = 0. The 
constant £^ — £^ can then be absorbed into a redefinition of gjj(iy ■ After these simplifications 
we are left with only two non-zero coefficients 

n /-sb n fsb 

§ C 7 = - 1 C ' = --^, (54) 

where we have defined Clr = CClr- 

The coupling of the Z' boson to quarks is not related to its coupling to leptons. There- 
fore tight constraints from semileptonic decays, as we encountered in the case of a modified 
Z coupling, can be avoided here by simply switching off the Z' coupling to leptons. Such 
"leptophobic" Z' bosons can for example appear in models with an Eq gauge symmetry (see 
e.g. ref. [§3]). Since leptophobic Z' bosons avoid detection via traditional Drell-Yan processes, 
their mass is much less constrained allowing for larger values of the parameter £. 

Besides constraints from hadronic B decays we have to face constraints from B s -B s mixing 
to which tree-level Z' exchange contributes. We find for the AB = 2-Hamiltonian: 



1 ^VRR _ 471-2 



d LR = ^f l^rl I^tI 7- (55) 
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In contrast to constraints from hadronic B decays, the B s -B s mixing constraint in the 
(Re(C| b ), Im(£| 6 )) -plane depends on the parameter £ determined by the coupling constant 
9u(i)' an d the Z' mass M//. It gets stronger for smaller £, i.e. for smaller gu(\y and larger 
Z' mass Mz'- This behaviour, which might seem counter-intuitive at first sight, has its origin 
in the dependence of the hadronic decays on the parameter combinations C,f b = £ C,f b . If one 
chooses smaller £ values, one needs larger values of the FCNC couplings (f b in order to obtain 
the same effects in the hadronic decays. Since the B s -B s mixing coefficients in (|55p depend 
quadratically on the (f b , this procedure sharpens their constraints. 



4.3.2 Results 

Considering eq. ([54"]) . one easily sees that the three scenarios C£ 6 7^ 0, (r an d C£ b = Cr ^ 
exactly correspond to the scenarios q? 7^ 0, q' 9 7^ and = q' 9 7^ in our model-independent 
analysis, except for a normalisation factor. In this way the exclusion regions from the 2 a 
constraints and the confidence levels from the fit can be immediately read off from figs. 0] 
and provided one rescales the axes by an appropriate normalisation factor and rotates the 
pictures by 180° to take into account the minus signs in eq. (|53p . 

In figs. [7] and [5] we present our results for the ££ 6 and the ( R b scenarios with the meanings 
of the green region and the red-hatched ring being the same as in the preceeding sections. In 
addition the 2a region for A s is shown for different values of £. We recognise that there is very 
little overlap of the region allowed by hadronic constraints with the region preferred by A s in 
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Br(B s ->*V) Br(B,-»p *) 




-2-1012 -2-1012 



~sb -sb 

Re(£ L )/k SM | Re(^ VksMl 



Figure 7: Enhancement factors of Br(B s 07r°) and Br(£? s —> 4>P°) for C£ h 7^ 0. The red-hatched ring 
corresponds to the SM uncertainty. The green area is allowed by the 2a constraints from all hadronic 
B decays while the regions inside the blue lines are compatible with the constraint from B s -B s mixing. 
From the biggest to the smallest region they stand for the cases £ = 1/10, £ = 1/25 and £ = 1/100, 
respectively. 



the (j£ case. The same holds for the ( L b = Cr scenario not shown. This behaviour is easily 
understood: The observables A^4cp an d A s both call for NP with a large imaginary part. The 
branching ratios of hadronic B decays depend linearly on the real part of Cl'r & t leading order, 
thus they pull the Cl'r values towards the imaginary axis. The observable A s , on the other 
hand, depends quadratically on Clr an< ^ favours values on the diagonal Re(C£ fe i j) = — lm(^ L b R ). 
For the Q L b setting this situation is relaxed due to the weak constraints from B — > ttK such 
that one can solve the two experimental discrepancies in A^cp an d in A s at the same time. 

From the diagrams we see further that the B s -B s mixing constraint is very tight. It prohibits 
large effects in B s — > cpiT ,(f>p for realistic values of the parameter £ < 1/25. For £ = 1/25, 
which would correspond for example to gu{iy ~ 9 and Mz' ~ 400GeV, and for £ = 1/50 we 
present the maximum enhancement factors in tab. [71 These numbers are obtained abandoning 
the la region of the B — > irK fit and requiring only agreement with the la constraints. We 
find that enhancement of a factor ~ 5 is possible in the ££ 6 and ^r scenarios whereas no effect 
can occur in the ( L b = ( R b case because of eq. ([33]) . For £ = 1/100 the constraints from B s -B s 
mixing become so strong that no effect in B s — > (j)iT , (ftp would be detectable. A measurement 
of a significant enhancement would therefore set a lower limit on £, equivalent to an upper limit 
on the Z' mass. 
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Br(B s ->*V) Br(B,-»p *) 




~sb -sb 

Re(^ B )/k SM | Re(^ fl VksMl 



Figure 8: Enhancement factors of Br(B s — > 4>ir°) and Br(S s — > for 7^ 0. The meaning of the 
coloured regions is the same as in fig. [7J 



4.4 MSSM 

Super symmetric effects in B decays have been studied in an enormous number of publications 
but most often hadronic decays have not been considered in such studies because of their large 
theoretical uncertainties. In the MSSM with conserved R-parity, all new flavour-changing in- 
teractions can be related to the squark mass matrices and enter all kinds of B decays via loops 
with virtual squarks and gauginos/higgsinos. Therefore one can expect supersymmetric contri- 
butions to be of roughly the same size in hadronic modes as they are in leptonic, semileptonic or 
radiative modes, so their effects will be most easily found where the uncertainties are smallest. 
However, once hints for supersymmetry are found in clean decay channels, one will also look for 
confirmations of these observations in other modes. Therefore we find it interesting to study 
the possible size of isospin-violation in the MSSM and whether large effects in the rare decays 
B s -> epp and B s -> 07T° can be expected or not, taking into account present experimental 
constraints on supersymmetric flavour- violation. Besides, we investigate whether the deviation 
of the AAcp measurement from the SM prediction can be explained in the MSSM, as it has 



been claimed recently [8l|, [8 



Throughout this chapter, we use the MSSM conventions of the SUSY Les Houches Accords 



(SLHA) 83j,[8j] and diagonalise the sparticle mass matrices exactly. We do not make use of the 
mass- insertion approximation (MIA), which means that we are not limited to the case where 
off-diagonal elements in the mass matrices are small with respect to the diagonal elements. 
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Scenario 


Br(B s ^7r u ) 


Br L (B 3 ^</>p u ) 


Br(B s -></>p u ) 


Br SM (B s -^7r°) 




Br SM (B s -></>p°) 








2.5 (4.9) 


3.6 (5.6) 


3.3 (4.9) 




i 

25 


^ 


3.6 (5.7) 


3.7 (5.7) 


3.4 (5.1) 






/■sb /-sb 

_ ^R 


1.0 (1.7) 


1.0 (1.7 ) 


1.1 (1.8) 








1.9 (3.7) 


2.6 (4.0) 


2.4 (3.6) 


e = 


1 

50 


Cf ^ o 


2.6 (4.1) 


2.6 (4.1) 


2.4 (3.7) 






/■sb /sb 
— Si? 


1.0 (1.7) 


1.0 (1.7) 


1.1 (1.8) 



Table 7: Maximal possible enhancement of the B s branching ratios compatible with the constraints 
from all other hadronic decays and B s -B s mixing for the Z' model. 



4.4.1 Flavour-violation in the down-squark sector 

First we consider the scenario where flavour-violation arises in the down-squark sector, i.e. 
where the 6x6 down-squark mass matrix contains flavour- violating elements. In this case we 
expect the largest SUSY effects in b — > s transitions to stem from gluino-(down-)squark loops 
since these loops come with the strong gauge coupling a s . Neutralino-squark loops arise from 
exactly the same off-diagonal elements but are suppressed by weak couplings, so they can be 
neglected to a first approximation. 



The authors of ref. [73J have analysed such a scenario and have found that a significant 
amount of isospin-violation can only occur via b — > sqq box diagrams with virtual gluinos and 
down squarks (q = u, d) and not from photon- or Z-penguin diagrams. A necessary condition 
therein is that the SUSY-breaking masses of the right-handed first-generation squarks, m~ and 



m~, are very different from each other. Recently, this idea has been seized in ref. [85J] and 
studied in the light of new data, especially measurements of CP asymmetries in B — > irK 
decays and of B s -B s oscillations. It was found that the B — > irK data can only be reproduced 
in a tiny region of the parameter space of the model. 



In contrast to these findings the authors of ref. [81[ have found a large impact of the 
gluino-mediated photon penguin in a mass-insertion calculation and state that sufficient isospin- 
violation is generated to explain the A^cp data. However, we find that this results from a 
missing factor — a/67r in eq. (41) of ref. [81[. We have performed a scan over the MSSM 
parameter space without using the mass-insertion approximation. Our full results for the 
AB = 1 Wilson coefficients are written down in Appendix (Dj We find the gluino-mediated 
photon penguin to yield corrections below the 3% level to the SM coefficients of the EW penguin 
operators for all the points which allow for a diagonalisation of the squark mass-matrices with 
eigenvalues greater than (100 GeV) 2 and satisfy the experimental constraints from b — > sj. 
Such corrections are negligible for the prediction of A^4cp- 

We conclude that we find no sizeable enhancement of EW penguins in the MSSM with 
flavour-violation in the down-sector. Neither can we explain the A^Icp discrepancy in this 
scenario nor can we expect large NP effects in the decays B s — > 0ir° and B s — > (j)p°. 
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4.4.2 Flavour-violation in the up-squark sector 

Supersymmetric flavour-violation can also arise in the up-sector via off-diagonal elements in 
the hermitian 6x6 up-squark mass matrix. In this scenario, penguin and box diagrams with 
virtual charginos and up-type squarks can provide sizeable contributions to B decays. We have 
calculated all of these diagrams, the results are given in Appendix iDl 

The small Yukawa couplings occuring in the quark-squark-chargino couplings and in the 
squark mass matrix strongly suppress the effect of certain off-diagonal elements of this matrix. 
The only relevant flavour-violating entries for b — > s transitions are therefore those correspond- 
ing to cx — £l and cl — tj? mixing. We define them via 


m| = m? [ 1 6& L * | 
,0 1 

and 

(M?) 62 = ^=m/) 32 = 6 u d L ml (57) 

with a generic squark mass mx. The remaining flavour-conserving elements are specified by 
m? = Diag(m?,m?,m? fl ), (Tjj)u = (^17)33^3 and tan/3 = 10. 

We have performed a scan over the free parameters in the value ranges defined in the table 
in fig. [10] to identify possible sources of large isospin- violation. In fig. [9] we plot the Z-penguin 
contribution to \qj\ and \q$\ {qi is defined in eq. (|3ip ) over the sum of all contributions to 
\qi\ and \qg\ for 30000 random points. We see that considering the Z-penguin only is a very 
good approximation for all points yielding large isospin-violation. In this way our scenario is 
essentially equivalent to the one with a left-handed flavour-changing Z coupling k s £ discussed in 
sec. !4.2l and we will stay in this approximation in the following. Since a non- vanishing k s £ breaks 



electroweak symmetry 2M ^ must involve the vacuum expectation values v Ut d- Therefore it 
is almost exclusively sensitive to £32^ and not to 832 L and we will neglect the latter in the 
following. 

Having calculated the MSSM mass spectrum for a given parameter point we apply the following 
constraints: 

• physical squark and chargino masses > 100 GeV , 

• Br(B — > Xgj) compatible with data at the 2a level, 

• chargino contribution \C^\ < |Cff| ~ 0.22. 

The last constraint ensures that fine-tuned points passing the Br(£? — > X s ^f) constraint are not 
considered. 



It has been suggested [82| that a non-vanishing can generate isospin-violating effects 
large enough to explain the A^lcp discrepancy within the QCDF framework. This would clearly 
be an interesting perspective for supersymmetric effects in purely isospin-violating decays. 
However, we cannot confirm this statement in our framework. In fig. [9] we find that Cg can 
be enhanced by about 25% and C7 by about 8% with respect to the SM in the presence of a 
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Figure 9: Left: (Right:) Z-pcnguin contribution to (|gg|) over the full {\q$\) for 30000 points in 
the MSSM parameter space. Blue: points excluded by the bound on |CV 7 |. Red: points excluded by 
Br(_B — > X s ^j). Orange: points passing all constraints. Not displayed: points excluded by lower bounds 
on squark and chargino masses. 



non- vanishing 6^ L Q. In section [3] we have seen that a 25% effect in Cg or an 8% effect in C-j 
are not enough to generate a large A^cp an d are also not sufficient to enhance the branching 
fractions of B s — > (fin and B s — > (fip° sizeably. A NP effect of this size would be hidden in the 
theoretical uncertainty and thus be unobservable. These findings are illustrated in figs. 1101 and 
[TTT where we display on the one hand a figure similar to fig. [3] in the complex k s £ plane and on 
the other hand a zoomed version of the upper plots in fig. [6] To both figures we add the k s £ 
values resulting from chargino-induced flavour-violating Z couplings in a scan over 30000 points 
in the MSSM parameter space as defined in the table in fig. [TO] to illustrate the statements of 
the last paragraph. 

We note that there is an important difference between our calculation and the one of [jil ]. 
namely the treatment of strong phases. We obtain all of these phases directly from QCDF, 
where they are suppressed by either a s (nib) or Aqcd /m B ■ Thus we typically find small strong 
phases even though their prediction comes with a large uncertainty. In the framework of 



ref. 82J only absolute values of penguin-to-tree ratios are predicted from QCDF whereas the 



corresponding strong phases can assume arbitrary values between and ir [871 ]. In this way, 
large CP asymmetries can be generated even without a large NP contribution and also the 
fact that A<jp(B~ — > ir°K~) and Acp(B° ~~ > tt + K~) have opposite sign is no longer puzzling 
because the phases of the various tree and penguin topologies are uncorrelated. In contrast, 
our calculation reduces the 2.5a SM discrepancy in AAqp only marginally. From fig. [TO] we 
can read off that the vast majority of the allowed points are outside the 2.2a region. Only a 
few rather fine-tuned points are between 2.2a and 2a. 



^Including the naturalness constraint |<532 iZ '| < 0.59 [§(| the effects become even smaller. 
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Figure 10: Discrepancy between A^4cp i n theory and experiment as a function of kI^/IksmI as defined in 
sec. 14.21 From light to dark the coloured regions denote 2.2a, 2a and la. On top we add the k^/I^SmI 
values resulting from chargino-induced flavour-violating Z couplings in a parameter scan as defined in 
the table. Blue (orange): Points (not) excluded by the bound on |CV 7 |. Not displayed: Points excluded 
by the lower bounds on SUSY masses. 



5 Conclusion 

In this article we have studied the possibility of probing isospin-violating NP in hadronic B 
decays. We have proposed to test the EW penguin sector of the effective weak Hamiltonian 
via the decays B s — > (j)7T° and B s — > c/)p and provided a detailed phenomenological analysis of 
these two modes in correlation to other hadronic B decays. 

Our analysis is motivated by discrepancies found in B — > irK decays, which are to date 
the best-measured hadronic b — > s decays. In particular, the 2.5 a discrepancy found in the 
observable AAcp can be interpreted as a sign of NP in the EW penguin sector of the theory. 
We have demonstrated in a model-independent analysis that this discrepancy can easily be 
resolved by an additional NP contribution to the EW penguin operators ,...,Q^q if it is 
of the same order of magnitude as the leading SM coefficient Cg M . An exception are parity- 
symmetric scenarios where the contributions to PP decays cancel. Whereas the solution in the 
case of NP in is, as expected, due to a new contribution to the colour-allowed EW penguin 
amplitude, we have pointed out that, in the case of NP in C*p , the solution mainly comes 
about via a weak annihilation contribution in the QCDF framework which has a surprisingly 
large imaginary part. In particular we have found for the case of equal new contributions to 
C-j and Cg that, even though these contributions tend to cancel in the EW penguin amplitude, 
the AAqp discrepancy can still be solved via the EW penguin annihilation amplitude, a fact 
that had not been noticed before. For various scenarios we have performed frequentist fits 
to B — > ttK data. We have found the fit to work well for NP in Cg while NP in Cy is 
only poorly constrained from B — > ttK alone. Especially in this case, the PV counterparts 
B — > pK and B — > ttK*, which carry a different interplay of chiralities, give valuable additional 
information. We have seen that present experimental data already set strong constraints: a 
new EW penguin amplitude can be at most ~ 5 times larger than the SM one; on the other 
hand a new amplitude of about the same size as the SM one is required to solve the AAqp 
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Figure 11: Zoom of the upper plots in fig. [5J On top we add the k£VI k sm| values resulting from 
chargino- induced flavour-violating Z couplings in a parameter scan as defined in the table in fig. 1101 with 
the same colour coding as in fig. 1101 



discrepancy. 

Whether the discrepancy in AAqp really is a manifestation of NP, however, is not clear 
because the experimental data are still not conclusive due to the large uncertainties in the 
theory prediction and the still somewhat low statistics of the measurements. The long-standing 
problem of large theoretical hadronic uncertainties in hadronic decays can, however, be partially 
addressed exploiting the large variety of hadronic decay channels. In the case at hand, the best 
test is provided by the B s — > (ftir and B s — > (ftp modes, which are purely isospin- violating 
decays and are dominated by the EW penguin amplitude. Since these decays are not related 
to other hadronic B decays via flavour symmetries, their analysis requires a determination of 
the hadronic matrix elements from first principles and we used the QCDF approach to this 
end. From the full QCDF results we have derived simple approximate expressions (eqs. (|2ip - 
(|'23|) . (|33|) ) which reproduce the B s — > (j)TT°, (ftp amplitudes with high accuracy for arbitrary 
scenarios with NP in Cj \ ...,c['q. By quoting these formulae we facilitate the study of these 
decays without an extensive implementation of the QCDF framework. 

In this work we have performed the first analysis of the impact of NP in EW penguins on 
B s —7- (ftir , (ftp . A new EW penguin amplitude of the same size as the SM one can easily enhance 
the B s —7- (ft-ir and B s — > (ftp branching ratios by an order of magnitude. We have performed 
a quantitative analysis parameterising NP in EW penguins in a model-independent way, at 
the level of Wilson coefficients, and studied the maximum enhancement of the B s — > (fin , (ftp 
branching ratios for various scenarios, with respect to the result of our B — > Kir fits and with 
respect to constraints from other hadronic B decays. The results displayed in tab. [3] show that 
in many cases a large enhancement is possible. Particular exceptions are parity-symmetric 
models which have no impact on the VV decay B s — > (ftp and scenarios with (approximately) 
equal contributions to and C^ 1 which cancel in B s — > (ftir° . 

A survey of concrete NP models has been performed in section [U where we have considered 
a modified Z° penguin, a model with an additional U(l)' gauge symmetry and the MSSM. 
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In such models additional constraints arise from the semileptonic decays B — > X s e + e~ and 
B — > K*l + l~ and from B s -B s mixing. In case of the modified-Z°-penguin scenario we have 
found that the semileptonic constraints still allow for NP to an extent which is sufficient to 
resolve the A^4cp discrepancy. On the other hand, they prevent the B s — > (J)ir°,(j)p decays 
from developing an enhancement which beats the hadronic uncertainties of the SM prediction. 
Therefore, a large effect measured in these decays would rule out the modified Z coupling. 
The semileptonic constraints can for example be avoided in a model with an additional Z' 
boson whose couplings to leptons can be switched off. Our analysis has shown that in this 
scenario constraints from hadronic B decays and B s -B s mixing can be fulfilled simultaneously 
only at the 2 a level. The tight constraints from B s -B s mixing limit a potential enhancement of 
B s —7- <^>vr° , (f>p° to a factor ~ 5. Furthermore the occurrence of such a measurable enhancement 
requires a large gu(iy coupling and/or a light Z' boson. Finally our conclusion for the MSSM 
is that it is impossible to obtain a new contribution larger than about 25% in Cg and about 
8% in CV, which is clearly not enough to generate a large A^4cp ° r a significant enhancement 
of the B s decays. The AAqp discrepancy can be reduced only marginally in this way. Only a 
few rather fine-tuned parameter points reduce it from 2.5a to 2a. 

We stress again that the decays B s — > (fiir , <pp° are highly sensitive to isospin- violating NP. 
Their measurement would complement the analysis of B — > Kn decays and could shed light 
on the "A^4cp puzzle" . For this reason we strongly encourage experimental efforts to measure 
these decays at LHCb and at future Super-B factories. 
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Appendices 

A General framework for the calculation of hadronic B decays 

Throughout this work, we study hadronic B decays into two light mesons using the framework 
of QCD factorisation U 17]. This framework is based on the well-known effective weak Hamil- 
tonian for AB = 1 transitions given in ([5]) and described in section lA.il Matrix elements of this 
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Hamiltonian are calculated in an expansion in Aqcd /m B ■ We give a few details concerning the 
generalisation of QCDF to a Hamiltonian containing operators with flipped chiralities in sec- 
tion |A2] and specify our input parameters in section lA.31 In section IA.4I we collect expressions 
for obtaining various observables from the QCDF-calculated amplitudes. 



A.l The effective weak Hamiltonian 



The starting point for the analysis of hadronic B decays is the parameterisation of high-energy 
transitions in terms of effective four-quark operators multiplied by short-distance coefficients. 
In the SM, this leads to the effective weak Hamiltonian in eq. ([5]). The short-distance coefficients 
Ci are calculated in the MS renormalisation scheme at the scale jUew ~ My/. Their low-scale 
values at /jew ~ n^b are obtained from renormalisation-group equations (RGE). In the SM, 
only C\ is 0(1), while C2 and the QCD penguin coefficients C3... 6 arise at order 0(a s ) and 
the electroweak penguin coefficients arise at order 0(a), albeit partly enhanced by factors of 
x t w = mt/ rn w an d/or 1/ sin 2 6w- The complete analytical expressions for these coefficients 
are written down e.g. in 66]. 



Besides the SM operators, we also consider the possibility that New Physics gives rise to 
the so-called "mirror" penguin operators Q[, obtained from Qi by a global exchange of left- 
and right-handed chiralities of the quark fields, (V — A) (V + A). We thus replace in eq. ([5]) 



(58) 



Given the enhancement of CV(mev 
modified RGE scheme presented in [IE 
as leading-order coefficients, i.e. 



) and Cg(^Ew) by xtw and/or 1/ sin 2 6w, we follow the 
for the SM coefficients and consider the enhanced terms 



C Q LO 



o 



(M W ) = - 



xtw 



+ 



3 sin 6w 



(\QB Q (x t w) - 4C Q (x t w)) 



(59) 



with the loop functions Bq(x) and Cq(x) from [66J. To be consistent we neglect at the same 
time electromagnetic corrections to the QCD penguin coefficients as well as any mixing of 
CV,...,Cio into Ci,...,Cq. Since this treatment improves the RGE evolution for C7, ...,Cio 
and since it is exactly these coefficients we are interested in, it is well suited for our analysis. 
By contrast, for the NP contributions we use the standard treatment for the leading-order 
RGE. 



A. 2 Operator matrix elements in QCDF 

From the effective Hamiltonian ([5|) the decay amplitude for the process B — > M\M<i can be 
calculated as 

A bIm iM2 = (M[ h h4 h) \n cS \B), (60) 

where h refers to a helicity amplitude in case of a decay into two vector mesons, B -> V1V2 
with h = 0, +, — . 

In QCDF the matrix elements of the effective Hamiltonian are organised in terms of flavour 
amplitudes on which are directly related to the topologies of the underlying weak transition, 
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for example colour-allowed EW penguin, colour-suppressed tree etc. Analogous amplitudes 
represent the corresponding weak annihilation transitions. The topological amplitudes on in 
turn can be decomposed into coefficients a, which are in direct correspondence to the operators 
Qi in the effective Hamiltonian. At NLO these coefficients a, contain the perturbative non- 
factorisable vertex-, penguin- and spectator-scattering corrections governed by the factorisation 
formula. For a complete description we refer to The expressions given there can easily be 
generalised to account for the mirror operators by adding a topological amplitude a[ to each 
of the a.i [S8| . A similar generalisation applies to the annihilation amplitudes . Here we only 
need to add the expressions for «3ew4EW which read 



a 



3EW 



(MiM 2 ) 



-a / /(MiM 2 ) + a / /(MiM 2 ), 
a'i{M l M 2 ) + a'/(MiM 2 ), 
(MiM 2 ) - Jf(M x M 2 ), 
(M 1 M 2 )-a' 7 p (M 1 M 2 ), 
(M 1 M 2 ) + a' 7 p (M 1 M 2 )) 



-'!> 



,'p 



if MiM 2 = PP, 
if MiM 2 = PV, 
if M 1 M 2 = VP, 
if MiM 2 = V°V°, 

if MiM 2 = y ± y ± , 



a 



'p 

4EW 



(MiM 2 ) 



'10 



ho 

ho 

-Jv 



(AfiAfj 
[M\M 2 \ 
{M\M 2 \ 



) - rf ^ di{M x M 2 
+ rf2 d£{M x M 2 ) 



x 

„M 2 
X 



^ a{(M 1 M 2 



I J± 



a? {M 1 M 2 ) + Tf> a'g{M 1 M 2 



,' p 
ho 



(MiM 2 ) + rf 2 4 P (M!M 2 )) , 



if MiM 2 
if M 1 M 2 
if MiM 2 
if M 1 M 2 
if MiM 2 



PP, 

py, 

VP, 

v°v°, 
v±v± 



(61) 



with f^ 1 = F^ Ml (0)/F£^ Ml (0) being a ratio of form factors, such that f^ 1 ~ m B /A QC D 
and Z^ 1 ~ Aqcd/^b in the heavy-quark limit 17]. The NLO expressions for the coefficients 
Ob i are equivalent to the ones for the coefficients a% from [16| up to the replacement C% — > C[. 
Additionally in the transverse amplitudes of B — > VV decays one has to flip the helicities, i.e. 
the expressions for af are needed for a[~ and vice versa. 

The pattern of signs appearing in eq. (|6ip is a consequence of the fact that matrix elements 
of the mirror operators are related to the SM ones by parity 



(M 1 M 2 \Q i \B) 



which implies that the amplitudes involve the coefficients a only in the combinations [5 



-ilMiM 2 (MiMalQ'jP), 
(') 



(62) 



o?(Mi,M 2 ) - VMl M 2 a i p (M 1 ,M 2 ). 



(63) 



Here r\M\M 2 = il is the parity of the final state. For PP and longitudinal VV final states, 
we have tjm 1 m 2 = 1 whereas for PV final states r}M x M 2 = ~ 1- I n this manner left-handed 
and right-handed NP give rise to different signatures and correlations among PP, PV or VV 
decays. Exploiting this feature can be very important in order to probe the chirality structure 
of a potential NP model. 



A. 3 Input parameters and tree to penguin ratios in B — > nK, B — y tttt 

In the framework of QCDF the decay amplitudes depend on quite a few input parameters such 
as form factors, Gegenbauer moments of light-cone distribution amplitudes, etc. In tab. [8] we 
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QCD scale and running quark masses [GeV] 



A- 



MS_ 



0.231 



m b (m b ) 



4.2 



m c (m b ) 



ffi.s 



1.3±0.2 0.090± 0.020 



m q /m s 



0.0413 



CKM parameters 



A 



0.225 



cb 



0.0415 ±0.0010 



\v ub /v c 



cb\ 



0.085; 



-0.025 
-0.015 



1 



(70 ± 10)° 



sin(2/3) 



0.673 ±0.23 



B meson parameters 



Lifetime 



Decay constant 



r[ps] 



/b [MeV] 
A b [MeV] 



B~ 



1.638 



1.525 



210 ± 20 



200 



+250 



1.472 



240 ± 20 



200 



+250 

-0 



Pseudoscalar-meson decay constants and Gegenbauer moments 



/[MeV] 
ai, ati t ± 



li 



131 


0.20 ±0.15 



K 



160 
0.06 ± 0.06 
0.20 ±0.15 



Vector-meson decay constants and Gegenbauer moments 



/[MeV] 
/-LpVfeV] 
"i,+ 

"2,02,+ 



209 ± 1 
165 ±9 


0.1 ±0.2 



218 ±4 
185 ± 10 
0.06 ±0.06 
0.1 ±0.2 



221 ±3 
186 ±9 


0±0.3 



Pseudoscalar-meson form factor at q = 



D 



7T 



0.25 ±0.05 



B ^ K 



0.34 ±0.05 



K 



0.31 ±0.05 



Vector-meson form factor at q = 



A 



B 



P 



o.3o; 



h0.07 
-0.03 

0.00 ±0.06 
0.55 ±0.06 



B^K* 



0.39 ±0.06 
0.00 ±0.06 
0.68 ± 0.07 



B K 



0.38185§ 
0.00 ± 0.06 



0.65 



+0.14 
-0.00 



Tabic 8: Summary of the theoretical input parameters for hadronic B decays into two light mesons. All 
scale-dependent quantities refer tou = 2 GeV unless otherwise stated. The values represent the most 
up-to-date values taken from [li, 17, 89j-93|. 
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provide a list of up-to-date values. We use the updated value of as obtained in |94j]; as 
of Ar = Ar we assume the lower value 200 MeV as suggested by exclusive hadronic decays, 
see 17, [HO, HI, Particular attention deserves the choice of Xh,a, i-e. the numbers which 
parameterise our ignorance of non-perturbative physics occurring in spectator scattering and 
weak annihilation due to the exchange of soft gluons. For Xh we use the definition of [la ]. 



m B 



X H = (l + p H e t +*)hi^- 



(64) 



with a non-perturbative scale A/j = 500 MeV. The default value is pu = and we estimate 
the uncertainty setting pn = 1 and freely varying (f>jj between and 2tt. In the light of 
experimental data, X\ requires more attention. The reasoning goes as follows. 





Figure 12: Experimental and theoretical values of the penguin-to-tree ratio defined in eq. (|65[) . For 
explanations see text. 
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A good test of the QCDF hypothesis is to look at the vrif-penguin to 7T7r-tree ratio 
which can be directly related to experimental observables: 



ai(7T7r) + a2(TT7r) 





V ub 


fn ' 






fl< . 



T{B- 



7T 



■K°) 



T(B- 



-7T°) 



1/2 



exp ' o ino +0 - 030 

u - iuu -0.018- 



(65) 



To a good approximation this ratio relates a pure QCD-penguin decay to a pure tree decay and 
allows to eliminate the uncertainty from the heavy-to-light form factor. Since the tree decay 
B~ —7- 7r~7r° suffers from small uncertainties in QCDF and is expected to receive negligible 
contributions from NP, eq. (|65p probes the accuracy of the QCDF prediction for the penguin 
amplitude a\(ixK\ In absence of sizeable NP contributions to QCD penguins it constrains 
the uncalculable weak-annihilation contribution which is responsible for the lion's share of 
the theoretical uncertainty in the penguin amplitude. Though subleading in Aqcd/^b, /3f is 
known to be numerically enhanced. 

Using the parameters in tab. [8] and expressing /3| via the complex 0(1) parameter Xa as 
ref. lfj with 

X A = {l + PA e i ^)hx r ^- (66) 



m 



with the default value pa = 0, we find 



ai^mr) + a2(7T7r) 



- U - U ' 8 -0.015> 



(67) 



which is slightly smaller than older results. Repeating the analysis of [16] for the PV and VV 
modes we obtain fig. [T2J The plots show the experimental central values and uncertainties of 
the ratio in eq. (|65p (circles around the origin) for B — > irK and three related decays. They 
are combined with "limagon" curves representing the corresponding theory predictions where 
the phase 4>a is freely varied between and 2ir while we set pA = 1 (1.5,2) for PP and PV 
modes and pa = 0.6 (1.0, 1.5) for VV modes to obtain the blue (purple, yellow) curves. The 
red dots correspond to pa = 0, 1, 1.5 (pa = 0, 0.6, 1.0 for VV modes) with (pA fixed as in the 
scenario "S4" in for PP and PV modes and (f) A = -40° for VV. 

These results lead us to the conclusion that, in the light of present data, we prefer to change 
the widely used treatment with pa = as default and pa = 1 (pa = 0.6 for VV) for the variation 
of (pA in order to have a more conservative estimate of the theoretical uncertainty. Nevertheless 
we confirm that pa < 2 is a reasonable upper bound for weak-annihilation contributions in 
QCDF. We adopt pa = 1-5 (pa = 1-0 for VV modes) as our default value, keeping the default 
for 4>a as above, and estimate the uncertainty with pa unchanged and 4>a G [0,2-7r). 



A. 4 Calculating observables in hadronic B decays 

Starting from a decay amplitude A(B — > M1M2) the corresponding decay rate can be calculated 

S 



as 



r(B MiM 2 ) 



167rm£ 



\A(B -> M X M 2 ) 



(68) 
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with a symmetry factor S. We have S = 1/2 if M\ and M2 are identical and S = 1 oth- 
erwise. For decays into two vector mesons, where three different helicity amplitudes ex- 
ist, one replaces \A(B -> M 1 M 2 )\ 2 by a sum over the three possible helicities of the fi- 
nal state, J2h=o - + \-A h (B — > V^l V2) | 2 . The branching ratios are easily calculated from 
Br(B ->■ MiM 2 ) = r(.B -> MiM 2 )/r to t with r to t being the total decay width of the B meson. 
Our theoretical predictions always refer to CP-averaged branching ratios defined as 

Bl(B -*T) = \ (Br(J5 -> /") + Br(£ /)) , (69) 

where /, / stand for the final state M\M<i and its CP-conjugated state. The CP-conjugated 
decay rate is calculated by replacing A in eq. ([68 p by the corresponding CP-conjugated am- 
plitude. This amounts to a change of sign for all the weak phases while strong phases are 
unchanged. 

The direct CP asymmetries read 

. Br(£? -> /) - Br(5 ->■ /) 

ACP = Br(S-/)+Br(5-/) (70) 

and the longitudinal polarisation fraction /l(B —> M1M2) is defined as 

, ,p , 11 M >_ l-^B^V^)! 2 .... 
/l( B -» M,M 2 ) - (71) 

B Isospin-violating observables in B — >■ 7r.Kr 

In this appendix we present the observables which are sensitive to isopin violation in the 
B — > nK decay modes and which we use to calculate our 2 a constraints. 



First, one has ratios of any two different decay rates. Using the parameterisation (jlOp and 
neglecting terms which are quadratic in the r{ as well as the annihilation contribution r^ w 
which has only a small real part, one has 6 different ratios which read 1(3, U] 



o Br(B- -> tt°K~) , P , 

R B C = 2 =J- ^J- -1 + 2 Re(r EW + rg w ) - 2 Re(r T + r c ) cos 7 , 

o 1 Br(B° -> vr+K") . r N 

R n = 2 B? ( go IT ^o) - l + 2Re(r EW + rg w )-2Re(r T + r c )cos 7 , 

R f = 2 ^- wh-o TlA ~ l + 2Re r EW -2Re r c cos 7 , 

r_ Br(i?° ->■ 

= ^H^^| * 1 + 2R ^T)co S7 - 2Re(,.C w ), 

Rl = ^ ^ B - *IJ -1-2 Re(r T + 2r ) cos 7+2 Re(2rE W + rg w ) , (72) 
t_ Br(£? u — > 7T U A U ) 
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where To and r_ are the lifetimes of the neutral and charged B mesons, respectively. NP in 
EW penguins as in (fTB"j) enters the ratios Rc,^ ,n through 

Re(rEw) -» Re(r E w) + Re(f E w) cos 5, 

Re(rg w ) -> Re(rg w ) + Re(r£ w ) cos 6. (73) 

Experimental data on the Rc,n^ ,7T can be used to constrain the NP contributions f E w an( A 
fg W . Note that the Rc,l^ ,n involve different combinations of f E w and r^ w an d thus they are 
sensitive to different linear combinations of the electroweak penguin coefficients ci , c['q . 
Therefore, it depends on the specific NP scenario in consideration which of the Rc,^'^ give the 
best constraints. 

Beyond being responsible for the universal QCD penguin contribution, isospin relations 
account for the approximate equation 

T(B~ -> tt~K°) - 2T(B~ -»• tt°K~) « 2T{B° -> ^°/?°) - r(5° -> ^+A~) (74) 



known as Lipkin sum rule [96J. In the strict isospin limit both sides of this equation vanish 
identically. This is reflected in the fact that R^ n in (|72p are both equal to one, apart from 
isospin-violating terms of order 0{ri). These linear terms are generated by the interference of 
the isospin-violating parts of the amplitude with the QCD penguin part. The special property 
of (|74p is now that these interference terms on the left- and righthand side of the approximate 
equation cancel each other. For this reason (f?4"|) can be used to construct a purely isospin- 
violating observable, namely 

r- Br(g° -> ^) ± rpIMB- - ,»*-) = 

r_Br(B° ^vr+K") + r r(B- ^ vr"K°) K%J y ' 

In a similar way it is possible to construct observables with a high sensitivity to isospin 
violation out of the direct CP asymmetries. To this end we consider the two differences 

A^p = A C p(B- -> tt°K-) - A C p(B -> tt + K-) = AA C p 

AA° CP = A CP {B- -> vr'A ) - Acp(£ -> vr A' ). (76) 



In the parameterisation (jlOp and in the presence of NP in electroweak penguins (JT3J) one 
finds AAqp = AA^p up to terms quadratic in the r,. The observable AA^,p is given in (|15p 
and represents the famous AAcp showing a 2.5<r discrepancy with current data. A precise 
measurement of AA^p could therfore shed light on the A^4cp discrepancy. 

Since 7r°K s is a CP-eigenstate into which both the B° and the B° meson can decay, we 
have mixing-induced CP violation in this decay channel. The corresponding observable Sqp is 
defined via the time-dependent CP asymmetry as 

where Ccp = — ^4cp is the direct CP asymmetry, up to a sign. Although S n x is not sensitive 
to isospin-violation in particular, it will be affected by a solution of the "A^4cp-puzzle" via a 
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NP contribution f^w- The reason is that tew has to come with a large new weak phase 5 in 
order to have substantial impact on Avlcp- Including the new electroweak contributions and 
neglecting Re(rg W ) we find 

S C p(B° -> 7T°K ) ~ sin 2/3 + 2Re (rc) cos 2/3 sin 7 — 2Re(fEw + ^ew) cos 2/3 sin 5. (78) 

We collect the theoretical and experimental results for the observables defined here, as well as 
the CP-averaged branching ratios and direct CP asymmetries in tab. [TJ 

Finally we give the QCDF expressions of the ratios used above: 



?'T 



?"EW 



x («) 



3A f 



a\{-KK) 
oc%{<kKY 



IKk 6c%(itK) 



rc 



a;; 



2 a%{nK) 



3 ^ C ,EW(^) 

2 d|(7riiC) 



(79) 



C The fit 

The basic idea of our fit is quite simple: we calculate within a NP scenario the expected values 
for a set of observables as a function of NP parameters qi. We then compare these values at 
each point of a grid in the {q{\ parameter space to experimental data. The points for which 
the experimental and the theoretical results are closest are most likely to be realised, i.e. they 
represent the (/j-values for which the theoretical prediction describes best the experimental 
measurements. Technically, this comparison is performed by evaluating at each grid point the 
X 2 function 

X 2 (tel) = E (Xjthe ° ({g 2 }) " XjeXp)2 ; (80) 

where the sum runs over different observables Xj. In this notation a^theo represents the theo- 
retical prediction of the observable and Xj eX p is the corresponding experimental mean value. 
<7jexp stands for the la experimental uncertainty (symmetric around the mean). 

The non-trivial part of the analysis is the implementation of the theoretical error. Here we 
follow the .Rfit scheme 0j. Our basic assumption is that experimental data approximatively 
yield a Gaussian distribution of an observable but a theoretical calculation does not. The 
latter depends on a set of input parameters like form factors, decay constants etc. for which no 
probability distribution is known. The Rht scheme corresponds to a frequentist approach and 
it assumes no particular distribution for the theory parameters, only that they are constrained 
to certain allowed ranges. The theoretical and experimental uncertainties are then combined 
in the following x 2 function: 

cxp— Xj thco|— CT j thco) 2 - f I _ I ^ 

-p. ^ l^jcxp -^j thco I Wjtheo; , . 

"jexp (81) 

otherwise. 
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Here we suppress the dependence on the parameters Since we often encounter asymmetric 
theory intervals, notated as (^theo)-^!^ 6 "';^ > we generalise eq. (f8Tj) to 



x 2 = E< 



hco,inf Zjcxp) -f - / _ \ 

32 11 J-j exp ^ V-^J thco u j theo, mi A 



(■^j theo ^Jtheo, inf ^jexp) 



(^j cxp ( x j thco~l~0"j thco, sup)) 

a 2 

3 cxp 



if Xj'exp > (^-j'theo "H thco, sup) j (^2) 



otherwise. 

We calculate <Tj theo, sup and cjj t heo, inf using the theory input given in tab. [8] and Appendix lA.31 
adding the resulting uncertainties in quadrature. 

Using the x 2 function in eq. (f82j) it is possible to define confidence levels (CLs) by means 
of the function 

CL({ 9i }) = ; / e-^t^/^dt, (83) 

v / 2^r(iV dof /2) 7a x2({(?! }) 

where Ax 2 is the x 2 function after subtraction of its minimum: A% 2 = X 2 — nhn(x 2 ). ^dof 
is the number of degrees of freedom of free model parameters. Setting CL = 1 — 68.27/100, 
CL = 1 - 95.45/100, CL = 1 - 99.73/100 and CL = 1 - 99.99/100 we find the 1 a, 2 a, 3 a and 
5 a confidence levels respectively. 

In our fits we only include quantities which are derived from independently measured ob- 
servables, i.e. AAq P , Ai^ p , ,Scp(-B — > ir K ) and one out of the three categories of ratios 
rB,k,tt ^ gee Appendix [B]) . Thereby we select the category which is most sensitive to the NP 
scenario under consideration. In this way we can on the one hand avoid to overweight a partic- 
ular observable in the fit and on the other hand avoid that the fit is pulled to large (^-values by 
discrepancies of quantities carrying only a small sensitivity on NP. We note here that, since the 
jjB,Kn are ra tj os f branching fractions, which we assume to be Gaussian in experiment, their 
probability distributions derived from experimental data are not exact Gaussians as required 
by eqs. (|81ll83p . Comparing fits to ratios of branching fractions to fits to the differences of 
the corresponding branching ratios (which follow a Gaussian distribution), we have checked 
that the qualitative outcome of the fits in terms of preferred regions in the complex gj-planes 
is not tarnished but that the contour lines are sharpened due to the reduction of theoretical 
uncertainties in the ratios. 



D SUSY contributions to penguin coefficients 

In this section we quote our results for the initial conditions of the Wilson coefficients C3 io, 
CV 7 and Cs g and of their mirror counterparts at the SUSY mass scale. We decompose the 
Wilson coefficients Cs,..,io into contributions C g from gluon penguins, C 7 from photon penguins 
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and C^ L1 , C^ L2 , C^ R1 , C^ R1 (q = u, d) from box diagrams: 



c 3 






= c 5 + 


1 ^tLL2 , /~<LL2 

2 U s„ > 


c 5 




c (i 


= c g + 


1^LR2 , ^<LR2 
2 B u ' Bd 


cv 


r< i r»LRi ^lri 
— o 7 -t- o Bu - u Bd , 


C s 


_ qLR2 


r<LR2 
~ U B d , 


Cg 


s"~i . /^yLLI ^LLl 

— u 7 + u Bu - o Sd , 


Cio 


_ qUJI 
B u 


~ <^B A i 



(84) 

The coefficients C" are obtained from these expressions by the replacements L f-> i? and 
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Cg T . In addition Z-penguins contribute to C\ according to eq. ([4Tj) . The individual 



contributions are obtained by calculating chargino-(up-)squark loops and gluino-(down-)squark 
loops. We do not consider neutralino-(down-)squark exchange since it is suppressed with re- 
spect to the gluino contributions involving the strong coupling. Our results are in agreement 
with similar calculations in the context of semileptonic B decays 97H99I] and with the gluino 
boxes in ref. 



We neglect additional operators arising from b —> sbb transitions because they contribute to 
B decays only in higher orders. Moreover we neglect box diagrams with more than one flavour- 
violating squark line. As for the squarks of the first two generations we assume approximate 
degeneracy of the corresponding elements of m^, likewise for m? and m~. Yukawa couplings 
of the first two generations are set to zero. In this way the box diagrams depend on common 
masses rrm L = mj^ { r mu R and ^j^) for the left-handed (right-handed) squarks of the first and 
second generation. 



We use the conventions of the SUSY Les Houches Accords 83j, |84| with only one exception: 
The CKM matrix is denoted as V whereas the chargino mixing matrices are named U and V 
instead of U and V. 



D.l Chargino contributions 

We use the loop functions written down in 98|, 99] and the mass ratio 

/ \ 2 

(85) 



\ Km / 

Generation indices i,j = 1,2,3, squark indices r,s = 1,...,6 and chargino indices m,n = 1,2 
are always summed over in the following. We abbreviate the chargino-quark-squark couplings 
by 

If im = W m2 KZ (5 ak + AU% ak ) V£ (% - AUl%) (<% + AC/*,) . (86) 

The matrices At/^, AC/]| and AC/^ account for chirally enhanced corrections to the tree-level 
quark-squark-chargino coupling. Explicit expressions for these matrices are given in Ref. (86[. 
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The quantity Y\ represents the modified Yukawa coupling incorporating tan /3-enhanced cor- 

rections (resummed to all orders in perturbation theory). An explicit formula for permit- 
ting complex SUSY mass parameters can be found in Ref. (lQOl ] . 

Neglecting Yukawa couplings of the first two generations we find that C3 10 = 0. The 
coefficients C^.^io are constructed from 

C l = Q/1 fo n T/ T/ * / ^ r r2»n r r3*m4 \ x u r c m ), (87) 

24V2GF7rV3 3 V32(m-+) z 

C<? = nz n 7} T/t ; ^r^mrrSm^i '(Vcm), (88) 

32V2GFvrl/33V 3 * 2 (m-+ ^ 

"^^ Am 

Ck L1 = — 1=^- m~+m-+U n iU*iT^2n^r^mDo(rn% ,mj. ,m 2 ~+ ,m 2 + ), (an) 



48a^V2G f .7rV r 33^2 ™ ' d '" ™ r ' Xm ' X " 

1 L L / 2 3 

V fc=l 

£ sr {2(7 (ml sJ m|+,m| + )m-+m 5 ,+W n iW^ 1 - C 2 (m? s ,m| + ,mJ+)V m iV* 1 } 

Am An -*. m Am Art J 



(91) 



The remaining box coefficients vanish. The magnetic and chromo-magnetic coefficients are 
given by 

= „ JkTi ta ta* / T^-^r2m ( ~T -^r3m^2 '(Vc™) ~~ ^3™^! ^(^M r c m ) J ; (92) 

^2G F V^Vz2\ m xV V m b J 

=—7=77— 7—7— r^r^m ( ^, mr r3*m^6 ( x «rC m ) ~ ^ r 3mK ( x u r c m ) ) ■ (93) 



D.2 Gluino contributions 

Here we use the loop functions written down in 9?], 98] and 

3 - 3x + (2 + x) log(x) 
Fio{x) = — 



These functions depend on the mass ratio 

Xgd a 



Writing 



G L si = KijfSji + AUtj,), 



12(x - If 



777 7 



G* = K d SJ+3 (8 ji + AU d R 



(94) 

(95) 
(96) 
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we split the 6x6 down-squark mixing matrix 1Z into a lefthanded and a righthanded 6x3 
block G L and G R . The matrices At/^ and AJTn account for chirally enhanced corrections to 
the tree-level quark-squark-gluino coupling [1011 ] . With a sum over all indices understood, we 
find (q = u, (ijil 

4\/2«a s F 6 (x„ d ), (97) 

a s 



qLR2 _ 



C v = o fnn tr t,* 2 G % G * ( C ^ 9 d s ) ~ C A F w (x gds )) , (98) 

^ = lOS^^VS {^l D o( m l m l m h<) + D2(mlmlml,ml)) , (99) 

= ~ 36^SS? 2 H^^' m i' m l'0 " 7^(mi,m|,m|. J mj t )) , (100) 

CbT = ~ 54 ^^ 2 BA)(mf,mf,mf r ,mfJ + 5£ 2 (mf,mf, m^O) , (101) 

E G ^G r tG^C 2 (m?,m| r ,m| s ). (103) 

i=l 

The magnetic and chromo-magnetic coefficients are 

- GjJGj,(C i ,F,(z sJ .) - (2C F - C^)F 2 ( %(i .))~) , 

with the group factors Cp = | and = 3. 

The corresponding expressions entering the Wilson coefficients C[ are obtained from 
eqs. (|97|) - (|105|) by the simple replacement L f-> with the exception of the Z-penguin which 
reads 

K « = -£ E G? 3 G^G£*G? 2 *C 2 (mJ, ro^m*). (106) 

i=l 



(104) 



(105) 



*The sign of D 2 in C^ L2 differs from ref. [73 
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